NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


DISSERTATION 

INTERACTION  OF  LASER  BEAMS  WITH 
RELATIVISTIC  ELECTRONS 

by 

Douglas  W.  Small 
March  1997 

Dissertation  Supervisor:  William  B.  Colson 

Approved  for  public  release;  distribution  is  unlimited. 

DTIC  QTJALTT^  INSPECTED  a 

19980223  123 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Publrc  ([^!^inqe^)^'*J^eX<inVthe  colleaion  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 

1.  AGENCY  USE  ONLY  (Leave  Hank)  2.  ^ORTDATE^  3.  ^OWTm^NI 

7  DATES  COVERED 

issertation 

4.  TITLE  AND  SUBTITLE 

INTERACTION  OF  LASER  BEAMS  WITH 

RELATIVISTIC  ELECTRONS 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Small,  Douglas  W. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943-5000 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect 
the  official  policy  or  position  of  the  Department  of  Defense  or  the  U.S. 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is 
unlimited . 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

Motivated  by  the  desire  to  put  a  free  electron  laser  (FEL)  weapon  on  a  ship, 
the  FEL  and  the  related  process  of  Compton  backscattering  are  studied.  The  theme 
of  the  majority  of  this  work  is  the  interaction  of  the  Gaussian  optical  mode  with 
a  beam  of  relativistic  electrons. 

Classical  FEL  theory  is  reviewed  in  Chapter  II.  Simulations  based  on  the 
classical  theory  are  used  in  Chapter  III  to  study  a  proposed  1  kW  (kilowatt) 
infrared  FEL.  In  Chapter  IV,  simulation  is  used  to  study  the  problem  of  electron 
beam/optical  mode  overlap  in  an  ultraviolet  (UV)  FEL.  A  new  concept,  the  FEL 
with  a  short  Rayleigh  length,  is  studied  in  Chapter  V.  The  idea  is  tested  on  the 

UV  FEL,  then  used  to  design  and  simulate  a  megawatt-class  FEL  for  ship  self- 
defense. 

An  analytical  calculation  of  the  Compton  backscattering  of  laser  light  is 
performed  in  Chapter  VI.  A  quantum  electrodynamics  (QED)  formalism  is  used  to 
find  the  spectrum  and  angular  distribution  of  photons  scattered  out  of  a  Gaussian 
optical  mode  by  relativistic  electrons. 

u  subject  terms 

free  electron  laser,  FEL,  simulation,  Compton  backscattering 

15.  NUMBER  OF  PAGES 

113 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

Standard  Form  298  (Rev  2-89) 


Prescribed  by  ANSI  Std  Z39-18 
298-102 


NSN  7540-01-280-5500 


i 


Author:  _ 

Approved  by: 


Approved  by: 
Approved  by: 


Approved  for  public  release;  distribution  is  unlimited. 

INTERACTION  OF  LASER  BEAMS  WITH 
RELATIVISTIC  ELECTRONS 

Douglas  W.  Small 
Lieutenant,  United  States  Navy 
B.S.,  Marquette  University,  1988 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

DOCTOR  OF  PHILOSOPHY  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1997 


'  Hai  inloo  \A /  Orvv 


Douglas  W.  Small 


William  B.  Colson 
Professor  of  Physics 
DissertatiOF  Supervisor. 


Robert  L.  Armstead  - 
Professor  of  Physics 


Karlheinz  E.  Woehler 
Professor  of  Physics 

Ron  J.  P 

ProfessoFf5f/9£dtriQfeVEngineering 


C=) — tL 

(/  James  H.  Luscoi 


Luscombe 
Professor  of  Physics 


1LX 


m 


Anthony/A.  Atchley/Chair,  Department  of  Physics 


Mi 


Haig  _L 


f.  P?E 

Ht 


A 


D.  Weir,  Associate  Provost  for  Instruction 


iii 


ABSTRACT 


Motivated  by  the  desire  to  put  a  free  electron  laser  (FEL)  weapon  on  a  ship, 
the  FEL  and  the  related  process  of  Compton  backscattering  are  studied.  The  theme  of 
the  majority  of  this  work  is  the  interaction  of  the  Gaussian  optical  mode  with  a  beam  of 
relativistic  electrons. 

Classical  FEL  theory  is  reviewed  in  Chapter  II.  Simulations  based  on  the 
classical  theory  are  used  in  Chapter  III  to  study  a  proposed  1  kW  (kilowatt)  infrared 
FEL.  In  Chapter  IV,  simulation  is  used  to  study  the  problem  of  electron  beam/optical 
mode  overlap  in  an  ultraviolet  (UV)  FEL.  A  new  concept,  the  FEL  with  a  short 
Rayleigh  length,  is  studied  in  Chapter  V.  The  idea  is  tested  on  the  UV  FEL,  then  used 
to  design  and  simulate  a  megawatt-class  FEL  for  ship  self-defense. 

An  analytical  calculation  of  the  Compton  backscattering  of  laser  light  is 
performed  in  Chapter  VI.  A  quantum  electrodynamics  (QED)  formalism  is  used  to  find 
the  spectrum  and  angular  distribution  of  photons  scattered  out  of  a  Gaussian  optical 
mode  by  relativistic  electrons. 
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I.  INTRODUCTION 


Today’s  US  Navy  is  the  most  powerful  the  world  has  ever  known.  It  can 
position  aircraft,  ships,  and  submarines  within  striking  distance  and  project  its  power 
deep  within  the  borders  of  almost  any  country.  Yet  today’s  ships  are  arguably  more 
vulnerable  to  attack  from  mines  and  anti-ship  cruise  missiles  (ASCMs)  than  ever 
before,  as  a  result  of  the  Navy’s  shift  from  blue  to  brown  water.  Ship  defense  against 
ASCMs  is  mainly  a  problem  of  reaction  time;  that  is,  there  is  not  much  time  (usually 
seconds,  depending  on  the  missile’s  speed  and  the  ship’s  detection  range)  between 
detection  of  the  missile  and  impact  of  the  missile  with  the  ship.  In  this  case  a  self- 
defense  weapon  which  travels  at  the  speed  of  light  would  be  most  useful.  With 
continuing  improvements  in  ASCM  speed,  stealth,  and  maneuverability,  the  problem 
gets  even  more  difficult  for  existing  defenses.  A  directed  energy  weapon  (DEW)  could 
be  the  solution. 

Experiments  in  missile  destruction  and  atmospheric  propagation  have  shown 
that  a  megawatt  (MW)  class  free  electron  laser  (FEL)  is  the  directed  energy  weapon  of 
choice  for  ship  self  defense  (SSD).  Megawatts  are  needed  in  order  to  deposit  the 
required  energy  on  the  target  in  a  short  amount  of  time.  More  power  probably  would 
not  shorten  this  time  because  thermal  blooming,  which  results  from  atmospheric 
absorption,  will  degrade  propagation.  The  FEL  is  needed  because  of  its  inherent 
tunability  and  designability,  characteristics  no  other  laser  system  can  boast.  The  laser 
weapon’s  wavelength  must  operate  in  very  narrow  atmospheric  absorption  windows  in 
order  to  propagate  the  beam  to  the  target,  avoiding  problems  like  thermal  blooming. 
The  FEL  has  another  advantage  over  other  laser  systems  in  that  its  gain  medium  is  a 
vacuum. 

While  SSD  is  seen  as  the  mission  area  where  laser  weapons  can  get  their 
start,  a  MW-class  FEL  aboard  a  ship  could  be  an  awesome  offensive  weapon  as  well. 
It  was  demonstrated  in  the  1970s  that  it  is  easier  for  the  laser  weapon  to  destroy 
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crossing  targets  than  incoming  targets.  It  is  possible  that  the  precision  and  power  of 
the  FEL  could  also  be  used  against  aircraft  and  their  pilots,  satellites,  ballistic  missiles, 
other  ships,  and  even  personnel.  Unlike  the  proposed  Arsenal  Ship,  which  merely 
packs  a  ship  with  more  missiles,  a  shipboard  directed  energy  weapon  would  truly  bring 
about  a  revolution  in  military  affairs,  which  would  hasten  the  obsolescence  of  the 
missile.  In  this  work  we  study  the  FEL  and  related  phenomena,  contributing  to  a  field 
of  scientific  research  with  the  goal  of  putting  a  DEW  on  a  ship. 

Several  topics  are  presented  on  the  interaction  of  laser  beams  with  relativistic 
electrons,  this  includes  the  free  electron  laser.  In  the  next  chapter,  an  overview  of 
classical  FEL  theory  necessary  for  understanding  the  ensuing  chapters  is  provided. 
We  describe  the  components  of  the  FEL  and  then  find  the  emission  wavelength  using 
basic  physical  arguments.  Using  the  Lorentz  force  and  Maxwell’s  equations,  we  derive 
a  set  of  two  equations  called  the  FEL  pendulum  and  wave  equations  which  describe 
the  FEL  interaction.  Using  these  equations,  we  find  expressions  and  describe  the 
mechanisms  for  FEL  gain.  Modifications  of  the  FEL  interaction  due  to  short  optical 
pulses  or  the  Gaussian  optical  mode  are  considered  in  the  last  two  sections  of  the 
chapter. 

Original  work  simulating  a  proposed  one  kilowatt  (IkW)  upgrade  to  the 
FIREFLY  FEL  at  Stanford  University  is  presented  in  Chapter  III.  This  upgrade  would 
represent  a  one  hundred-fold  increase  in  FEL  average  power  over  that  which  currently 
exists.  This  work  was  originally  presented  at  the  16th  International  FEL  Conference 
and  published  as  "Simulations  of  the  Stanford  FIREFLY  IkW  free  electron  laser"  in 
Nuclear  Instruments  and  Methods  in  Physics  Research  (NIM)  A  358,  ABS  44,  (1995). 

In  Chapter  IV  is  presented  original  research  concerning  optical  mode 
transverse  effects  in  ultraviolet  free  electron  lasers  (UVFEL).  Specifically,  it  is  the 
UVFEL  proposed  by  the  Laser  Processing  Consortium  (LPC),  a  group  based  at  the 
Thomas  Jefferson  National  Accelerator  Facility  (TJNAF),  that  is  simulated.  The  results 
were  presented  at  the  17th  International  FEL  Conference  and  published  as 
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"Transverse  effects  in  UV  FELs"  in  NIM  A  375,  ABS  61,  (1996). 

New  and  original  work  on  a  new  concept  in  FEL  design,  the  short  Rayleigh 
length  FEL,  is  presented  in  Chapter  V.  The  concept  is  first  tested  on  the  LPC  UVFEL. 
Those  results  were  presented  at  the  18th  International  FEL  Conference  and  will  be 
published  in  1 997.  We  then  use  this  concept  to  design  and  simulate  a  MW-class  FEL 
oscillator  suitable  for  ship  self-defense.  This  design  and  the  simulation  results  were 
presented  at  a  workshop  at  TJNAF  entitled  "Navy  MW-Class  SSD  FEL  Concepts" 
which  was  sponsored  by  the  Navy  High  Energy  Laser  office  at  the  Space  and  Naval 
Warfare  Command  (SPAWAR). 

In  Chapter  VI  we  present  new  and  original  calculations  of  Compton 
backscattering  of  laser  light  from  electrons.  This  process  has  been  proposed  as  a 
means  of  generating  x-rays  for  many  applications,  including  medical  imaging,  and  y- 
rays  for  a  so-called  y-y  collider.  Our  calculations  are  the  first  attempt  to  incorporate  the 
finite  size  of  the  laser  beam  into  the  usual  scattering  problem. 
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II.  CLASSICAL  FREE  ELECTRON  LASER  THEORY 


A.  THE  FREE  ELECTRON  LASER 

The  Free  Electron  Laser  is  a  device  which  transforms  the  kinetic  energy  of  a 
beam  of  relativistic  electrons  into  coherent  radiation  (Ref.  1).  A  transverse,  static 
magnetic  field  causes  the  electrons  to  wiggle  and  thereby  radiate  electromagnetic 
energy  (Ref.  1).  The  radiation  can  either  be  stored  in  an  optical  cavity  and  amplified 
over  many  passes,  as  in  the  FEL  oscillator,  or  amplified  in  a  single  pass,  as  in  the  FEL 
amplifier  (Ref.  1).  In  this  work,  we  will  be  concerned  only  with  the  oscillator 
configuration. 

The  major  components  of  the  FEL  oscillator  are  the  electron  accelerator, 
undulator,  and  optics.  The  accelerator  delivers  a  beam  of  relativistic  electrons  to  the 
FEL  with  energies  anywhere  from  a  few  MeV  to  GeV  depending  on  the  wavelength  of 
laser  light  desired.  There  are  many  different  types  of  accelerators  suitable  for  use  in 
an  FEL,  but  we  will  only  be  concerned  with  the  radio-frequency  linear  accelerator  (RF 
LINAC).  The  RF  LINAC  accelerates  pulses  of  electrons,  called  micropulses,  using 
either  standing  or  travelling  waves  of  an  electromagnetic  field  (Ref.  2).  One 
micropulse  can  be  accelerated  within  each  wavelength  of  the  RF  field,  called  a 
"bucket.”  For  example,  if  the  frequency  of  the  RF  radiation  is  f  =  500  MHz,  buckets 
are  separated  by  a  distance  c/f  =60  cm.  Electron  micropulses  are  created  and 
inserted  into  the  accelerator  by  the  injector  (or  electron  gun).  In  our  example,  if  the 
injector  also  operated  at  500  MHz,  it  would  fill  every  RF  bucket  and  micropulses  would 
be  separated  by  60  cm.  The  current  in  the  micropulse  is  called  the  peak  current  /, 
and  the  average  current  is  calculated  using  the  duty  factor  (the  ratio  of  the  micropulse 
length  to  the  distance  between  buckets).  If  we  take  /  =  500  A  and  a  micropulse  length 
le  =  0.6  mm  (le/c  =  2  ps).  then  the  duty  factor  in  our  example  accelerator  is  0.001. 
The  average  current  is  then  7  =  500x0.001  A  =  0.5  A. 
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The  undulator  is  a  set  of  permanent  magnets  or  electromagnets  arranged  to 
provide  a  periodic,  transverse,  magnetic  field  of  wavelength  X0  along  which  the 
electrons  traverse.  The  polarization  of  the  undulator’s  magnetic  field  can  be  either 
helical  or  linear.  The  world’s  first  FEL  used  a  helical  undulator,  but  the  linear 
undulator  is  more  common  now  (Ref.  3).  Undulators  can  range  in  length  from  1  m  to 
25m. 

The  mirrors  of  the  optical  cavity  (in  the  oscillator)  store  the  light  which  is  then 
amplified  over  successive  passes  of  the  electrons  through  the  undulator.  Total  optical 
cavity  losses,  including  the  outcoupled  light,  are  quantified  by  the  quality  factor  Q,  so 
that  losses  per  round  trip  of  the  optical  pulse  are  (100/Q)%  (Ref.  3).  The  optical  cavity 
is  generally  centered  on  and  several  times  longer  than  the  undulator. 

B.  THE  RESONANCE  CONDITION 

With  a  basic  knowledge  of  the  physics  of  the  FEL  interaction  we  can  define 
some  important  parameters  and  find  the  wavelength  of  the  emitted  radiation. 

Depicted  in  Figure  1  is  the  FEL  interaction:  the  interaction  of  relativistic 
electrons  with  the  undulator  magnetic  field,  But  and  the  radiated  electric,  £?r,  and 
magnetic,  Br,  fields.  At  the  top  is  a  schematic  of  an  FEL  showing  the  electron  beam 
path  through  the  undulator,  the  orientation  of  the  undulator  magnets,  and  the  mirrors. 
The  middle  sketch  is  a  blow-up  of  one  wavelength  \0  of  the  undulator.  The  undulator 
magnetic  field  Bu  oscillates  in  the  yz-plane,  causing  the  electrons  to  wiggle  in  the  xz- 
plane.  At  the  bottom  of  the  figure  is  a  blow-up  of  one  half  of  one  wavelength  X  of  the 
radiation  field.  A  portion  of  the  electron  trajectory  is  drawn  at  an  angle  to  the  z-axis.  It 
is  the  transverse  component  of  the  electron  velocity,  that  which  projects  on  the  xy- 
plane,  which  couples  to  the  radiation  electric  field,  allowing  energy  to  be  transferred 
between  them. 
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Figure  1 .  Relativistic  electrons  interacting  with  an 
undulator  magnetic  field,  Bu,  and  radiated  electric,  Er, 
and  magnetic,  Br,  fields  (after  Ref.  4). 


In  order  for  the  radiation  in  the  oscillator  cavity  to  grow,  the  electrons  must  lose 
energy  to  the  light  wave  on  each  pass  through  the  undulator.  The  Lorentz  force 
shows  that  this  energy  loss  is  dy/dt  $  •  %  where  $=V/c,  V  is  the  particle  velocity 
and  the  relativistic  factor  y  =  (1  -  f^2)-172  (Ref.  5).  Since  £:r  is  transverse,  that  is,  in  a 
plane  perpendicular  to  the  z-axis,  only  a  transverse  component  of  ^  couples  to  the 
cavity  radiation.  Driven  by  the  undulator  magnetic  field,  that  transverse  component  is 
of  the  form  pi  «  R e(Kelk°z)  where  K  is  a  constant  called  the  "undulator  parameter," 
ko  =  2ti/X0,  and  X0  is  the  undulator  wavelength.  The  electric  field  is  of  the  form 
Er  oc  Re(E0e,(te_b)f+<l>))  where  k  =  2k/X,  X  is  the  radiation  wavelength,  to  =  kc,  and  <f>  is 
an  optical  phase  factor.  Therefore,  we  find  that  dy/dt  <*  Re(KE0e'l{k+k°)z  ~  It  is 

evident  from  these  preliminary  considerations  that  the  electron  "ponderomotive"  phase 
defined  as 

C  =  (k  +  k0)z  -  cot  (1) 

is  fundamental  to  the  FEL  interaction.  For  relativistic  FELs,  this  phase  represents  the 
position  of  the  electron  within  a  radiation  wavelength. 

The  time  It  takes  the  electrons  to  traverse  the  undulator  is  Lic$z  =  lo-8. 
Rather  than  work  with  time  scales  on  the  order  of  10-8  seconds,  we  introduce  a 
dimensionless  time  t,  defined  as 


where  L  is  the  undulator  length  and  cpz  =  vz  =  c  in  the  relativistic  limit.  As  the 
electrons  traverse  the  undulator  x  ranges  from  0  to  1.  We  can  now  define  the  electron 
phase  velocity  as 

^)=^-  =  LHk  +  k0)^~k]  .  (3) 

The  phase  velocity  measures  the  rate  of  change  of  the  electron’s  phase  relationship 
(Q  to  the  optical  field. 
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In  order  that  the  rate  of  transfer  of  energy  from  the  electrons  to  the  radiation 
field  is  maximized,  the  wiggling  electron  and  the  electric  field  must  be  nearly  resonant. 
Looking  at  Figure  1 ,  we  can  explain  this  as  follows:  when  the  electric  field  is  pointed 
in  the  x-direction,  the  electron  trajectory  needs  to  have  a  component  in  the  x-direction. 
When  the  electric  field  phase  increases  by  n  and  is  pointed  in  the  (-x)-direction,  the 
electron  trajectory  must  also  have  a  component  in  the  (-x)-direction.  This  will  be  the 
case  if  exactly  one  wavelength  of  light  passes  over  the  electron  every  time  the 
electron  traverses  exactly  one  wavelength  of  the  undulator.  This  is  the  "resonance 
condition."  If  one  wavelength  of  light  passes  over  the  electron  for  every  undulator 
period  the  electron  traverses,  the  electron  will  maintain  the  same  phase  relationship 
with  the  light.  Resonance  thus  implies  that  the  phase  velocity  v(x),  must  be  zero. 


Invoking  the  resonance  condition,  we  set  v(t)  =  0  in  (3)  and  solve  for  X  with  the 
result  X  =  X0(1  -  pz)/pz.  From  the  definition  of  y,  we  find  y-2  =  1  -  p2  -  pf.  In  the 
helical  undulator  (as  determined  later),  p2  =  K2/ y2.  Solving  for  pz  we  get 


Pz  = 1  -  1r#  + 
2 y2 


(4) 


for  y  »  1 .  Substituting  (4)  into  the  expression  for  X  above  we  find 


(5> 
2  r 

which  is  the  wavelength  of  light  emitted  by  the  electrons  in  the  FEL.  This  expression 
shows  that  we  can  change  the  wavelength  by  changing  the  electron  beam  energy 
y me2  or  the  undulator  characteristics  Xq  or  K. 

In  the  proceeding  discussion  we  defined  and  then  discussed  resonance  for  a 
single  electron.  In  a  real  electron  beam,  there  are  about  1010  electrons  in  the 
micropulse;  all  of  them  could  not  possibly  be  exactly  resonant.  Although  we  define 

o  o 

resonance  as  v  =  £  =  0,  where  (  )  represents  differentiation  with  respect  to  i,  we  will 

find  that  in  the  low  gain  FEL  (section  E)  it  is  desirable  to  operate  at  just  above 
resonance  v(0)  =  v0  =  2.6.  We  can  understand  this  by  considering  that  about  106 
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electrons  are  uniformly  spread  over  the  length  of  an  optical  wavelength.  If  they  are 
resonant,  they  maintain  the  same  phase  relationship  with  the  light  throughout  the 
interaction  even  though  the  light  is  passing  over  them.  Thus  at  resonance  the 
radiation  from  the  ith  electron  at  £,■  will  exactly  cancel  that  from  an  electron  at  a  phase 
(C,  -  k).  This  cancellation  will  occur  for  each  electron.  To  avoid  this  cancellation,  for 
net  radiation,  the  electrons  must  radiate  together,  or  "bunch"  on  the  scale  of  the 

.  o  oo 

optical  wavelength.  This  requires  that  £  >  0  and  £  >  0.  We  will  discuss  this 
microscopic  bunching  in  part  three  of  the  next  section. 

C.  ELECTRON  DYNAMICS 

Electrons  in  the  undulator  interact  with  both  the  undulator  field  and  the  radiation 
field  which  has  been  set  up  in  the  cavity.  In  subsection  1,  we  find  the  electron 
trajectories  in  the  undulator  using  the  Lorentz  force  equation  of  motion.  We  then  use 
these  trajectories  in  the  second  subsection  to  deduce  measures  of  electron  beam 
quality.  In  the  third  subsection  we  assume  that  electron  radiation  has  set  up  a 
classical  radiation  field  in  the  optical  cavity  and  then  allow  it  to  interact  with  the 
electrons.  This  will  result  in  the  pendulum  equation,  which  describes  bunching. 

1.  Electron  Trajectories  in  the  Undulator 

A  linear  undulator  polarized  in  the  y-direction  causes  electrons  to  wiggle  in  the 
xz-plane.  In  addition,  since  electron  beams  have  a  small  radius  and  angular  spread 
(which  depends  on  the  quality  of  the  beam  from  the  accelerator),  electrons  off  axis  in 
the  y-direction  are  in  stronger  fields  and  are  therefore  focussed  back  toward  the  axis. 
This  results  in  so-called  betatron  oscillations. 

The  magnetic  field  of  a  linear  undulator  with  no  external  focussing  in  the  x- 
direction  is  given  by 

S  =  [Bx,  By  B2)  =  B0[ 0,  sinfkozjcoshf/foy),  cos(/c0z)sinh(fr0y)]  (6) 
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where  B0  is  the  amplitude  of  the  undulator  field  (Ref.  3).  Consider  an  idealized 
electron  which  is  injected  into  the  undulator  perfectly  (i.e.,  x(0)  =  px(0)  =  0)  in  the  x- 
direction.  The  electron  motion  is  governed  by  the  relativistic  Lorentz  force  equation  of 
motion  in  the  undulator  field  (Ref.  5): 


(7) 


where  e  =  lei  is  the  magnitude  of  the  electron  charge,  and  the  electron’s  momentum 
is  p  -  ymc$,  m  is  the  electron  mass,  and  c  is  the  speed  of  light  in  vacuum.  Note  that 
cgs  units  are  used.  Substituting  (6)  into  (7),  noting  that  dy/dt  =  0  since  the  magnetic 
field  can  do  no  work  on  the  electron,  we  get  the  set  of  equations 


eB0 

Px  =  :^[Min(/foz(0)cosh(k0y(f))  -  pycos(/c0z(f))sinh(/c0y(f))]  (8a) 

ymc 

P y  =  ■^Pxcos(/c0z(f))sinh(k0y(f))  (8b) 

fiz  =  -^PxSin(M(f))cosh(/c0y(f))  (8c) 


where  (  '  )  =  d/dt.  The  right  side  of  (8a)  is  an  exact  time  derivative  and  can  be 
integrated  immediately  with  the  result: 


Px  =  — —cos(k0z)cosh(k0y)  (9a) 

where  K  =  eBXQ/2nmcz  and  B  =  Bq/^2.  is  the  root-mean-square  (rms)  magnetic  field 
strength.  The  dimensionless  parameter  K  is  the  undulator  parameter  and  is  an 
important  quantity  in  FEL  physics,  as  shown  in  (5).  The  constant  of  integration  has 
been  set  to  zero  indicating  our  assumption  of  perfect  injection  in  this  direction.  The 
time  dependence  of  the  dynamical  variables  has  been  suppressed  for  brevity. 
Substituting  (9a)  into  (8b)  and  (8c)  we  get 


cK2k 

— j-2-cos2(/c0z)sinh(2/c0y) 
T 


(9b) 
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cK2kr 


-sin(2/c0z)cosh2(/c0y)  . 


Equations  (9)  are  coupled,  nonlinear  and  too  difficult  to  solve  analytically.  In 
order  to  investigate  the  xz-plane  motion  we  will  now  add  the  restriction  of  perfect 
injection  in  the  y-direction  also.  That  leaves  us  with 


Px  =  ~~~~cos(k  qZ  ) 


(10a) 


Py(0)  =  y(0)  =  o 

■  cK2k0 
Pz  =  — ~s\n{2k0z) 
T 


(10b) 


(10c) 


Now  in  the  FEL,  fC~1  and  ys>  1,  so  that  Kly  1  (Ref.  3).  This  allows  us  to  expand 
equations  (10)  for  the  trajectories,  keeping  only  terms  to  order  {Kly)2.  To  zeroth  order 
in  {Kly),  z{t)  =  c$zt  where  cj5z  we  take  to  be  the  average  electron  velocity  in  an 
undulator  period.  Using  (4),  we  find  cpz  =  c(1  -  (1  +K2)/2f  +  ...)  =  c.  Substituting 
z(t)  =  c$2t  into  (10c)  to  find  the  next  higher  order  contribution  we  get 

z(0  =  cM  - -^sin(2k0cM)  +  (11) 

4/Cq't 

which  is  already  to  order  {Kly)2.  To  find  x{t),  we  will  substitute  the  first  term  of  (11) 
into  (10a)  since  using  the  second  would  lead  to  higher  order  terms  in  x{t).  We  get 

yfOtS 

*(0  =  --^sin(MPz0+  (12) 

where  in  both  (11)  and  (12)  we  have  set  =  1  in  the  denominator. 

We  can  rewrite  (11)  and  (12)  in  terms  of  the  dimensionless  time  t  (see  (2))  as 


,  v  KX o 

x(x)  =  — ~~s\n{2KNx)  +  •  • 
<2ny 

z(t)  =  P7Lt  -  — ^sin(4n/VT)  + 


(13a) 


(13b) 
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where  N  =  L/XQ  is  the  number  of  undulator  periods.  As  expected,  the  x-component  of 
the  trajectory  is  the  wiggling  motion  caused  by  the  undulator  with  one  oscillation  per 
undulator  period.  The  z-component  contains  the  average  velocity  term  and  a  term 
with  a  fast  oscillation  (twice  per  undulator  period).  The  oscillation  is  superimposed  on 
the  average  motion  because  as  an  electron  with  constant  speed  makes  an  excursion 
to  the  side,  the  z-component  of  its  velocity  must  slow,  leading  to  an  oscillation  in  z- 
trajectory.  An  electron  makes  this  excursion  twice  per  linear  undulator  period,  which 
causes  the  z-trajectory  to  oscillate  twice  per  undulator  period. 

Now  that  we  have  found  the  trajectories  in  the  x-  and  z-directions  (13),  we  will 
find  the  motion  in  the  y-direction.  The  y-trajectories  are  called  betatron  oscillations 
and  are  caused  by  the  focussing  action  of  the  undulator  magnetic  field.  We  need  to 
use  a  more  realistic  beam  in  this  direction,  but  maintain  a  small  enough  beam  such 
that  k0y«:  1  and  sinh(2k0y)  =  2 k0y.  If  we  then  re-write  (9b)  in  terms  of  x,  the  result  is 


oo  2  K2L2k£  „ 

y  (x)  = - 5 — — cos2(2jrA/x)y(x) 

r 


(14) 


Since  a  typical  betatron  oscillation  takes  several  undulator  periods  to  complete,  we  can 
average  (14)  over  several  periods  (Ref.  3).  We  are  then  left  with 


°°/x  K2Lzk$  /%  aVp(y) 
y(t)= — — 


(15) 


which  is  the  equation  of  motion  for  a  simple  harmonic  oscillator.  The  solution,  and 
thus  the  y-trajectory,  is  sinusoidal.  The  harmonic  potential  is  l/p(y)  =  topy2^  where 
cop  =  KLk^y  =  2nNK/y  is  the  dimensionless  betatron  frequency.  In  a  typical  undulator 
N  =  100,  K  =  1.  y  =  100,  and  we  get  cop  =  2k,  or  one  betatron  oscillation  along  the 
undulator. 


We  have  only  treated  betatron  motion  in  one  direction,  but  there  can  be 
betatron  focussing  in  the  other  direction  as  well.  This  could  be  provided  by  external 
focussing  magnets,  or  by  machining  parabolic  pole-faces  on  the  undulator  magnets 


13 


(Ref.  3).  Neither  the  amount  of  focussing  nor  the  electron  beam  quality  have  to  be 
the  same  in  both  directions. 


2.  Electron  Beam  Quality  and  Emittance 

By  inspection  of  (15)  and  knowledge  of  Hamiltonian  mechanics,  we  know  the 
harmonic  betatron  oscillations  have  a  constant  of  motion 
HP  =  [®py2^)  +  y2(x)]/2  =  [copy2(0)  +  y2(0)]/2.  It  can  also  be  shown  that  when  an 
electron  is  injected  off  axis  (y(0)  =  y0)  or  with  an  initial  angle  (y(0)  =  LQy  where  Qy  is 
the  angle  from  the  z-axis),  its  z-velocity  is  decreased  by  Apz  =  -Hp/L2  relative  to  a 
perfectly  injected  electron  (Ref.  3).  Taking  differentials  of  (3)  we  find  Av  =  L/rApz  since 
k  s>  kQ.  Therefore,  an  imperfectly  injected  electron  experiences  a  decrease  in  phase 
velocity  given  by 

Av  =  (K2k$y  o2  +  ^e2)  (1 6) 

where  we  have  used  (5). 

In  a  real  electron  beam  from  an  RF  LINAC  there  can  be  =1010  electrons  per 
micropulse,  each  of  which  has  a  phase  velocity  change  away  from  the  average  given 
by  (16).  Numerically  we  don’t  simulate  all  1010  electrons,  but  instead  distribute  many 
electrons  in  both  injection  position  and  angle.  They  are  distributed  as  a  Gaussian, 
which  has  the  distribution  function 


fjq)-  exp(-q2/2a2) 
>/2n  a 


(17) 


where  a  is  the  standard  deviation  and  q  is  the  variable  being  distributed  (Ref.  3). 
Using  (16)  as  a  guide  we  define  the  standard  deviation  of  the  angular  distribution  as 


o«  = 


-  ZnA/y2^2 


(18) 


1  +  K2 

where  6  is  the  root  mean  square  (rms)  initial  angle  from  the  z-axis.  The  standard 
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deviation  of  the  position  distribution,  also  called  the  dimensionless  electron  beam 
radius,  is 


o*=Virr«  <i9> 

where  T0  is  the  rms  initial  distance  from  the  z-axis  (Ref.  3).  The  significance  of  the 
factor  VXD i,  the  characteristic  optical  mode  size,  will  be  discussed  in  section  G  of  this 
chapter.  From  (1 6)  we  see  that 


®ve 


2%NKZI<qTq 
1  +  K2 


(20) 


is  the  standard  deviation  of  an  exponential  spread  in  phase  velocity  caused  by  the 
initial  Gaussian  spread  in  position  characterized  by  (19).  The  Gaussian  spread  in 
injection  angles  characterized  by  (18)  also  leads  to  an  exponential  distribution  in 
phase  velocity  with  the  same  standard  deviation  (Ref.  3).  The  exponential  distribution 
arises  because  the  phase  velocity  of  an  electron  decreases  no  matter  what  the 
direction  of  the  (Gaussian  distributed)  initial  position  or  angle.  The  distribution  function 
for  the  exponential  is  given  by 

f(q)  =  ,<o  (21) 


where  q  is  the  variable  being  distributed,  i.e.,  v/t  the  phase  velocity  of  the  Ith  electron 
(Ref.  3). 

The  injection  angles  and  positions  of  the  electrons  in  the  beam  are  governed 
by  a  quantity  called  the  "emittance,"  which  is  a  measure  of  the  quality  of  the  electron 
beam.  The  rms  emittance  is  defined  as  =  rcro0,  and  is  a  fixed  quantity  throughout 
the  entire  length  of  the  undulator  (Ref.  2).  That  is,  one  cannot  Change  the  radius  of 
the  electron  beam  without  affecting  its  angular  spread  and  vice  versa.  Solving  (18) 
and  (19)  for  e  and  T0  and  substituting  these  into  the  definition  of  emittance  yields 

=  VX  (22) 
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where  we  have  defined  the  dimensionless  emittance  e  as 

e  =  aeVoe  ■  (23) 

This  dimensionless  emittance  will  be  used  extensively  in  Chapters  IV  and  V. 


3.  The  Pendulum  Equation 

We  now  consider  the  interaction  of  the  electrons  with  both  the  undulator  and 
radiation  fields.  We  have  already  found  the  macroscopic  electron  trajectories  in  the 
undulator.  Our  goal  in  this  subsection  is  to  find  an  expression  which  describes  the 
evolution  of  bunching.  We  will  use  the  field  of  a  helical  undulator  in  this  section  as  it 
makes  the  problem  mathematically  easier.  At  the  end  we  will  state  a  prescription  for 
transforming  the  pendulum  equation  to  that  for  a  linear  undulator. 

The  magnetic  field  on  the  axis  of  a  helical  undulator  is  given  by 

3,  =  B0[cos(k0z),  sin (k0z),  0]  (24) 

and  the  radiation  field  present  in  the  cavity  can  be  represented  by 

K  =  E(x,y,z,t)[cosy,  -siny,  0]  (25a) 

&r  =  E(x,y,z,t)[sm\\r,  cosy,  0]  (25b) 

where  E(x,  y,  z,  t)  is  the  wave  amplitude  and  y  =  kz  -  wf  +  <f>  (Ref.  3). 

To  describe  this  problem  of  the  interaction  of  electrons  with  the  undulator  and 
radiation  fields,  there  are  five  coupled  equations: 


-f  =  -e(e  +  j!x0) 


(26) 


*1  =  .  £ 
dt  mcp 


y- 


vm 


(27) 

(28) 


but  only  four  unknowns:  Px,  Py,  P*,  and  y,  so  we  can  neglect  the  z-component  of  (26). 
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The  fields  £:  and  E$  in  (26)  and  (27)  are  the  sum  of  the  radiation  and  undulator  fields. 
Substituting  (24)  and  equations  (25)  into  the  x-  and  y-components  of  (26),  we  get 


-^( VSi)  =  E(  1  -  pz)[cos\)/,  -sin y,  0]  +  pzS0[-sin(/r0^),  cos (k0z),  0] 


(29) 


where  we  have  used  pi  =  ymc$x  and  E  =  E(x,y,z,t).  The  right  side  of  (29)  contains 
two  terms,  the  first  of  which  is  proportional  to  (1  -  pz)  and  is  surely  small  for  relativistic 
electrons.  We  will  neglect  this  term  and  keep  only  the  dominant  undulator  term  (i.e., 
that  which  is  proportional  to  B0).  We  can  then  integrate  (29)  and  set  the  constant  of 
integration  to  zero  indicating  perfect  injection.  The  result  is 


=  -j[cos(k0z),  sin {kQz),  0]  (30) 

where  K  =  eB^ZTimc2  is  the  undulator  parameter  since  B  =  B0  in  the  helical 
undulator. 

As  mentioned  earlier,  Z:  is  transverse  and  couples  with  only  the  transverse 
components  of  jl  that  is,  (30).  Therefore  substituting  (30)  into  (27)  we  are  left  with 

Y=^cos(C  +  <t>)  (31) 

y  me 

where  £  is  given  by  (1)  and  <t>  is  the  phase  of  the  optical  field  as  discussed  above. 

As  we  noted  earlier,  (28)  can  be  arranged  as  y~2  =  1  -  P f  -  Pf.  If  we 
substitute  (30)  for  px2  we  get 


(32) 


which  shows  that  a  change  in  y  does  not  result  in  a  change  in  pA,  only  in  pz.  Taking 
the  time  denvative  of  both  sides  of  (32)  results  in 


y  tMl 

1  1  +  K2 


(33) 


Construction  of  y/yfrom  (31)  and  equating  that  with  (33)  results  in 
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(34) 


r  _  1  +  K2  (k  +  k0)c  eEK  1% 

- fc—y5£c"«:+*) 

where  we  have  used  pz  =  tj{k  +  k0)c.  We  will  now  make  approximations  using 
k  »  k0  and  pz  =  1 ,  and  substitute  for  k  using  (5)  to  get 


••  2k0eEK 

C  =  — p - cos(C  +  <>) 

'fm 


(35) 


We  can  make  (35)  dimensionless  by  setting  d2/dt2  =  (c/L)2d2/dx2,  which  leaves  us 
with 


C  =  la  lcos(C  +  <|))  (36) 

which  is  the  pendulum  equation.  We  have  defined  the  dimensionless  optical  field 
amplitude  la  I  as 


la  |  _  4nNeKLE 
■fmc2 


(37) 


The  pendulum  equation  appropriate  for  the  linear  undulator  can  be  found  by  the 
substitution  K  -*  K[J0(Q  -  J&)]  in  (37)  (Ref.  3).  Here  J0(5)  and  J&)  are  Bessel 
functions  of  the  first  kind  of  order  0  and  1  and  £  =  K2/ 2(1  +  K2). 

oo 

Equation  (36)  shows  that  electrons  with  phase  -n/2  <  (£  +  ty)  <  nJ2  have  L,  >  0 

OO 

while  those  with  n/2  <  (£  +  <|>)  <  2nl2  have  £  <  0.  This  difference  in  phase 
acceleration  shows  how  electrons  which  are  initially  distributed  randomly  from 
(^  +  <j>)  =  -jt/2  to  (£  +  4>)  =  3ji/2  are  driven  together,  i.e.,  bunched,  near  (£  +  <)>)  =  jt/2  by 
interaction.  In  the  next  section,  we  will  see  what  affect  these  bunched 
electrons  have  on  the  optical  field. 


D.  THE  FEL  OPTICAL  WAVE  EQUATION 

Knowing  through  the  pendulum  equation  that  the  electrons  will  be  bunched  by 
the  FEL  interaction,  we  now  find  what  affect  those  bunched  electrons  have  on  the 
classical  radiation  field.  We  use  Maxwell’s  equations  driven  by  the  electron  current  to 
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determine  the  form  of  the  FEL  optical  wave  equation  (Ref.  3). 

Electrons  in  the  undulator  emit  spontaneous  radiation  because  they  are  made 
to  accelerate.  After  approximately  one  pass,  there  is  enough  radiation  in  the  optical 
cavity  to  take  on  a  classical  form  (Ref.  3).  Over  the  next  several  passes  this 
spontaneous  emission  spectrum  is  narrowed.  In  many  cases,  the  final  FEL  emission 
spectrum,  or  linewidth,  is  determined  by  the  length  of  its  optical  pulses  through  the 
Fourier  transform  (Ref.  3).  Short  optical  pulses  result  in  broader  spectra  than  long 
optical  pulses.  At  this  point  in  the  operation  of  the  FEL,  the  radiation  wavelength  is 
given  by  (5).  Since  we  know  that  the  emission  spectrum  of  the  FEL  is  narrow  (it  is  a 
coherent  light  source),  it  must  be  that  the  amplitude  and  phase  of  the  radiation  field  in 
the  cavity  are  slowly  varying  in  both  time  and  space.  We  will  invoke  the  slowly-varying 
amplitude  and  phase  approximation  in  the  following  derivation  of  the  FEL  optical  wave 
equation  (Ref.  3). 

The  fields  in  the  cavity  of  an  FEL  with  a  helical  undulator  are  given  by 
equations  (25).  Since  the  amplitude  and  phase  are  slowly  varying  in  time,  we  can 
write  the  vector  potential  for  these  fields  as 


where  t  =  [-/,  1,0]  is  the  polarization  vector  and  \y  =  kz  -  cot  +  <J>.  The  potential  for 
the  physical  fields  is  given  by  the  real  part  of  (38).  The  inhomogeneous  wave 
equation  for  /f  is  given  by 


(39) 


where  J  is  the  driving  current  density  (Ref.  5).  The  slowly-varying  amplitude  and 
phase  approximation  can  be  written  as  inequalities:  E  «  kE,  <>'  «  kp,  E  «  toE  and 
♦  «  oxf>  with  (  )’  =  d(  )/'6z  (Ref.  3).  In  the  wave  equation,  this  approximation  has  the 
effect  of  favoring  first  derivatives  over  all  other  terms.  Substituting  (38)  into  (39)  and 
carrying  out  the  z-  and  t-  derivatives  we  are  left  with 
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(40) 


Vf  +  2 ik^~  +  E(x,y,z,t)e^x’y'z'V  =  _il ™J.&e-i(kz  -  at) 

OZ  C  ut 

where  =  d2/dx2  +  d2/dy2  and  6*  is  the  complex  conjugate  of  e. 

The  driving  current  i/  is  given  by 

•A* ’V’z)  =  ~ec  X  8(x  -  */)  8(y  -  y,)  8(z  -  z,)  ft  =  -ec  2  8(3)(f  - ft)  ft-  (41 ) 

i  i 

where  8(x)  is  the  Dirac  delta  function,  T  is  the  position  vector  of  the  field  point 
(•*>  y i  z)>  7/  is  the  position  vector  of  the  i^1  electron,  and  the  index  /  runs  over  all 
electrons  in  the  micropulse.  Substituting  (30)  and  (41)  into  (40)  and  dividing  by  2 ik 
results  in 


l-~vf  +  !■#-  + 

4  n  + 


dz  c  dt 


Ee‘*  =  -2neK  £  - —  S(3)(f  -  Tj) 
i  Y i 


where  we  have  used  (1)  and  E  =  E(x,y,z,t). 


We  made  the  left  side  of  (42)  slowly  varying  by  invoking  the  slowly  varying 
amplitude  and  phase  approximation.  The  right  side  is  not  slowly  varying:  it  contains 
many  Fourier  components  due  to  the  presence  of  the  delta  functions  which  are  each 
infinitely  broadband.  This  means  that  while  the  left  side  describes  an  optical  field  that 
varies  slowly,  the  right  side  can  describe  an  optical  field  which  varies  arbitrarily  fast. 
To  smooth  the  right  side  we  must  replace  the  sum  over  all  electrons  by  an  average 
over  some  sample  electrons.  Since  the  fields  in  the  cavity  are  slowly  varying  over 
several  wavelengths  by  assumption,  we  only  need  to  sample  the  electrons  in  one 
section  of  the  micropulse  one  wavelength  of  light  long  and  use  the  local  charge 
density  to  weight  the  average.  The  result  is 


__L-V 


~~jT,  Ee"  =  “2 neK  p(z  -  c$zt)  < 


where  p(z  -  cp,f)  is  the  charge  density  of  the  traveling  electron  pulse,  and  <  •••  >  is 
an  average  over  sampled  electrons. 
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We  can  write  (43)  in  a  much  simpler  form  if  we  change  variables  to  2  =  z  -  ct 
and  use  the  dimensionless  time  x.  The  coordinate  2  then  follows  a  point  on  the 
optical  pulse  traveling  at  speed  c  from  which  the  electron  pulse  slips  back  at  speed 
c(1  -  (3Z).  In  this  coordinate  system  (let  z  ->  2)  we  find  that  L (d/dz  +  Vcd/dt)  =  d/dx. 
We  can  make  (43)  dimensionless  by  making  the  replacement  x,y  Vjr/XT  x,y  where 
'I’kL/n  is  the  characteristic  optical  mode  size.  The  result  is 


dx 


Ee'+  =  -2ti eKL  p  <  —  > 
Y 


(44) 


where  we  have  assumed  a  long  electron  pulse  so  that  p  can  be  considered  constant. 
Also,  if  the  FEL  efficiency  is  small  so  that  changes  in  y  are  small,  (44)  takes  the  form 


-V?  + 


_a_ 

dx 


a  =  -j<e 


-'Tv 


(45) 


where  a  =  la  le'<t>,  la  I  is  given  by  (37),  and  the  dimensionless  current  density  j  is 


8A/(e7tKI)2p 

-fmc2 


(46) 


To  get  the  wave  equation  for  an  FEL  with  a  linear  undulator,  we  make  the  same 
substitution  as  before,  K  -4  K(J0(fy  -  J^)),  in  both  (37)  and  (46)  (Ref.  3). 


Now  we  will  consider  the  physics  contained  in  (45).  First,  we  see  that  the 
growth  of  the  field  amplitude  91a  l/9x  depends  on  both  the  current  j  and  the  wavefront 
curvature  through  the  diffraction  operator  Vj.  This  unnecessarily  complicates  our 
discussion,  so  we  will  assume  for  the  moment  that  the  light  is  constructed  of  plane 
waves  so  that  E(x,y,z,t)  —>  E(z,t)  and  Vfa  -*  0.  Then  (45)  takes  the  simple  form 

a  =  .  (47) 

o 

We  see  that  if  j  =  0,  then  a  =  0  and  the  field  does  not  grow.  Also,  if  there  is  no 
bunching,  the  electrons  are  randomly  distributed  throughout  each  wavelength  of  light, 
<exp (-/£)>  =  0  and  again  the  field  does  not  grow.  Equating  the  real  and  imaginary 
parts  of  (47)  leads  to 
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31a  I  .  x. 

-fo-  =  -y<cos(^  +  4>)> 


(48a) 


|a|f  =  y<sin(C^)>  (48b) 

which  shows  that  the  field  amplitude  la  I  grows  when  the  electrons  are  bunched  near 
(C,  +  <)))  =  k  whereas  the  phase  <|>  is  driven  when  bunching  is  near  (£  +  <)>)  =  nl 2. 

The  wave  equation  (45)  and  the  pendulum  equation  (36),  when  coupled,  allow 
the  field  and  the  electron  phase  space  positions  to  evolve  self-consistently.  The 
pendulum  equation  describes  the  bunching  of  electrons  in  the  presence  of  the 
undulator  and  radiation  fields,  while  the  wave  equation  shows  how  much  the  radiating 
bunch  changes  the  optical  field.  The  two  equations  form  a  feedback  loop  which  under 
the  right  starting  conditions  leads  to  gain,  our  next  topic. 


E.  GAIN 


Now  that  we  know  how  the  electrons  affect  the  optical  field  and  vice  versa  we 
can  translate  this  interaction  to  gain.  FELs  can  be  separated  into  two  regions:  low 
current/low  gain  defined  by  j  1  and  high  current/high  gain  defined  by  j  »  1  (Ref.  3). 
As  it  turns  out,  the  low  gain  behavior  we  will  analyze  has  a  broader  range  of 
applicability,  up  to  at  least  j-  1.  An  FEL  with  y  =  100  will  exhibit  high  gain 
characteristics.  The  difference  between  the  low  and  high  gain  regions  is  in  the  growth 
rate  of  the  optical  field  a  (see  (47)).  We  will  find  that  the  gain  mechanisms  in  these 
two  regions  are  different  and  therefore  warrant  a  distinction  between  the  two. 

Gain  for  either  high  or  low  current  is  defined  as 


Lost  by  Electrons 


Optical  Field  Energy 


(49) 


where  the  last  equality  arises  from  conservation  of  energy.  Here  a0  is  the  optical  field 
strength  seen  by  the  electrons  upon  entering  the  undulator.  Using  (3),  (4),  and  (5)  it  is 
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not  difficult  to  show  that  the  change  in  electron  phase  velocity  due  to  a  change  in 
energy  Ay/y  or  wavelength  emitted  AAA  is 


Av  =  4nN— 
Y 


(50a) 


or 


Av  =  2kN 


AX 


(50b) 


so  that  the  energy  lost  by  one  electron  can  be  expressed  as 


AE  = 


(vp  ~  v) 
4%N 


ymc2 


(51) 


where  E  =  ymc 2  is  the  electron  energy.  The  number  of  electrons  in  a  volume  of  the 
optical  mode  dV  is  p FdV,  where  the  filling  factor  F  is  the  ratio  of  the  cross  sectional 
areas  of  the  electron  beam  and  optical  mode.  The  optical  energy  contained  in  that 
same  volume  dV  is  E2dV/4n  (Ref.  5).  After  some  algebra  we  can  write  (49),  letting 
j  ->  jF,  as 


G  = 


~^"<v0— v> 

a0 


(52) 


which  is  valid  for  both  low  and  high  current.  An  average  over  the  phase  velocities  v( 
in  the  sampled  portion  of  the  electron  beam  is  indicated  by  <  —  >. 

In  the  two  subsections  which  follow  we  will  discuss  gain  in  the  low  current  and 
high  current  FELs,  respectively.  We  will  find  expressions  for  the  single  pass  gain  in 
both  cases  as  a  function  of  x.  In  FEL  design,  low  gain  usually  implies  an  oscillator 
configuration  while  high  gain  (e.g.  j  =  104)  implies  an  amplifier  in  which  electrons 
execute  only  one  pass  before  the  light  is  outcoupled.  In  this  case,  the  dimensionless 
power  P  =  la!2  evolves  as  P  =  (1  +  G)P0  and  a  large  G  ensures  P  »  P0.  After  n 
passes  in  the  oscillator  Pn  =  (1  +  G)"P0  so  although  G  is  small,  over  many  passes 
the  light  is  amplified  until  P„  »  P0. 
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In  the  following  analysis  the  wave  equation  is  in  the  form  (47).  In  section  G  we 
will  discuss  the  effects  of  diffraction  on  the  FEL  interaction. 

1 .  The  Low  Gain  FEL 

In  exact  analogy  with  a  simple  pendulum,  the  tool  used  to  analyze  the  slow 
electron  dynamics  (bunching)  in  the  FEL  is  their  (C,v)  phase  space,  one  2n  section  of 
which  is  one  wavelength  of  light  long. 

In  the  low  gain  FEL,  j  «  1  implies  through  (47)  that  changes  in  the  field  are 
small.  Multiplying  both  sides  of  (36)  by  £  and  rearranging  terms  leads  to 

r  ■> 

la  lsin(C  +  4»)  =0  (53) 

and  we  recall  that  v  =  C,.  This  has  a  simple  solution 

r  *i  1/2 

v(Q  =  ±  2Hq  +  2\ a  lsin(C  +  <t>)|  (54) 

where  H0  =  v0/2  -  la  lsin(£0  +  (j>)  is  a  constant  of  the  motion  for  each  electron  starting 
at  (Co.  vo)  and  determines  which  phase  space  path  v(Q  it  will  traverse.  These  paths 
are  plotted  in  Figure  2.  For  electrons  starting  on  resonance  v0  =  0,  (54)  predicts  fixed 
points  at  (Co  +  <t>)  =  M2  and  3n/2.  Analysis  shows  that  the  fixed  point  at  k/2  is  stable 
and  surrounded  nearby  in  phase  space  by  harmonic  orbits  (Ref.  6).  The  fixed  point  at 

3ji/2  is  unstable  and  defines  the  path  characterized  by  H0  =  la  I,  called  the  separatrix, 
which  is 

vs  =  ±\'2lal  1  +  sin(C  +  4>)  1/2  .  (55) 

The  separatrix  separates  phase  space  into  regions  of  closed  and  open  orbits  and  has 
a  peak  to  peak  height  of  4 v  ia  I  (Ref.  6).  It  is  plotted  in  bold  in  Figure  2. 

In  Figure  3  we  have  distributed  20  resonant  (i.e.,  v0  =  0)  electrons  uniformly 
through  phase  space  and  show  their  evolution  as  x  :  0  -»  1 .  At  x  =  0,  the  electron 
position  is  gray  and  steadily  darkens  as  x  — »  1 .  The  phase  space  simulation  solves  the 


wave  and  pendulum  equations  numerically  with  j  =  1  and  initial  field  amplitude  a0  =  3. 


FEL  PHASE  SPACE  PATHS 
a=l  <J)=0 


Figure  2.  The  phase  space  paths  v(Q  traversed  by  the 
electrons.  The  separatrix  is  drawn  bold. 


Since  the  electrons  are  resonant,  they  start  and  remain  inside  the  separatrix  on  closed 
orbits.  One  can  clearly  see  in  Figure  3  that  the  electrons  become  bunched  at  £  =  rc/2, 
as  predicted  at  the  end  of  section  C.  Even  though  the  electrons  are  bunched, 
however,  gain  is  nearly  zero  because  the  electrons  are  bunched  at  the  wrong  phase 
for  field  amplitude  growth.  We  can  see  in  the  figure  that  as  many  electrons  move  up 
in  phase  space  and  gain  energy  from  the  light  wave  (see  (50a))  as  move  down  in 
phase  space  and  lose  energy  to  the  light  wave.  Something  must  be  done  to  upset  this 
balance  in  order  to  get  a  net  gain. 

It  is  interesting  to  note  that  Figure  3  shows  a  very  small  amount  of  gain.  This 
is  due  to  the  slight  phase  shift  given  by  (48b)  which  causes  the  bunch  to  move  just 
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slightly  toward  (C,  +  <(>)  =  n. 


Figure  3.  Phase  space  simulation  of  low  gain  FEL  with 
electrons  starting  on  resonance.  Gain  is  negligible. 


To  get  an  analytic  solution  for  the  gain  in  a  low  current  FEL,  expand  (36)  in 
weak  fields  la  I  =  a0  «:  k  with  the  result 

v(t)  =  v0  +  sin(£0  +  v0x)  -  sin(£0)J  +  (56) 

a0?  f  i  r 

lcos^2^0  +  2vot)  ~  cos(2Co)J  +  cos(v0t)  -  1  -  v0TSin(£0)cos(^0  +  v0x)  +  •  •  •  , 

where  a0  =  fa(0)'  and  0(0)  =  0  (Ref.  3).  When  we  perform  the  average  in  (52)  only 
the  term  proportional  to  a$  in  (56)  survives  and  the  result  is 

.  [  “  2cos(v0t)  -  v0xsin(v0x) 

G(T)  =  y  - - -  .  (57) 

v0 
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We  see  that  the  gain  at  the  end  of  the  undulator  (t=  1)  depends  only  on  the  initial 
phase  velocity  and  current  density  when  the  system  is  in  weak  optical  fields. 

In  Figure  4  is  the  result  of  a  simulation  which  solves  the  pendulum  and  wave 
equations  numerically  for  many  different  initial  phase  velocities  and  plots  the  final  gain 
G(  1=  1)  and  phase  shift  <|)(t  =  1)  versus  v0.  The  v0  axis  can  be  viewed  as  either  a 
function  of  electron  energy  or  emitted  wavelength  through  equations  (50).  In  the 
figure,  v0  ranges  from  -12  to  12,  j  =  1,  and  a0  =  0.001  ck.  The  gain  spectrum  is 
anti-symmetric  in  agreement  with  (57)  and  has  a  maximum  of  G  =  0.13  at  v0  =  2.6. 
The  phase  curve  shows  that  the  shift  is  maximum  for  resonant  electrons,  as  discussed 
above,  and  is  quite  small  where  the  gain  is  large. 


****  Gain  and  Phase  Curves  **** 
j=l  ao=0„001 

0.13 

0.0 

-0 . 13 
0.08 

0.0 

_ _ _ _ ! _ _ _ _ _ 1-0.08 

-12  v0  12 

Figure  4.  Gain  and  phase  spectra  for  low  gain  FEL. 

Gain  is  maximum  at  v0  =  2.6. 
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We  have  in  Figure  5  another  phase  space  simulation  using  the  same 
parameters  as  in  Figure  3  except  that  v0  =  2.6  in  order  to  maximize  gain.  In  this  case 
only  about  half  of  the  electrons  are  trapped  in  the  separatrix  while  the  other  half  are 
situated  on  open  orbits.  A  clear  bunch  of  electrons  can  be  seen  at  £  =  7t  which  affords 
maximum  gain  with  a  small  phase  shift  <|>.  Compared  with  Figure  3,  more  electrons 
lost  energy  to  the  light  wave  than  absorbed  energy  from  the  light  wave.  It  is  also  true 
that  if  the  electrons  started  at  v0  =  -2.6  the  result  would  be  maximum  absorption  as 
more  electrons  would  be  absorbing  energy  than  losing  energy.  This  behavior  could 
have  been  predicted  from  studying  Figure  2:  since  all  electrons  travel  nearly  the  same 
distance  along  their  orbit  (from  (36)),  starting  them  higher  in  phase  space  (i.e.,  v0  >  0) 
results  in  electrons  with  £  71/2  traveling  more  or  less  parallel  to  the  £  axis,  while 
those  with  >  n/2  travel  down,  losing  energy. 


***  FEL  Phase 

Space 

Evolution  *** 

3=1 

V3 

vo=2.6 

0.13 


-Tt/2 


Figure  5.  Phase  space  simulation  of  low  gain  FEL  with 
electrons  starting  at  v0  =  2.6  for  nearly  maximum  gain  of 
G  =  0.12  with  a0  =  3. 
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2.  The  High  Gain  FEL 

In  the  high  gain  FEL,  j  »  1  and  changes  in  the  optical  field  are  not  small. 
Energy  conservation  cannot  be  used  to  determine  the  gain  in  this  case  because  there 
is  a  substantial  change  in  both  the  optical  field  amplitude  and  phase.  It  is  still  possible 
to  obtain  an  analytic  expression  for  the  gain  in  weak  optical  fields;  we  merely  state  the 
result: 

G(x)  =  •^-exp[(-~-y')1/3V3x]  (58) 

where  v0  is  set  to  zero  (Ref.  3).  Unlike  the  low  gain  FEL  (57),  here  the  gain  is 
exponential  in  x.  We  can  ascertain  the  dependence  on  v0  by  using  the  same  gain  and 
phase  simulation  used  in  Figure  4.  Shown  in  Figure  6  is  that  simulation  with  j  =  100 
and  a0  =  0.001  <s  %. 


-12  Vo  12 

Figure  6.  Gain  and  phase  spectra  for  high  gain  FEL. 
Gain  is  maximum  at  v0  =  1 .4,  but  is  still  substantial  at 
v0  =  0  where  the  phase  shift  <(>  =  k/2. 
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The  maximum  gain  is  G  =  71  at  v0  =  1.4,  closer  to  resonance  and  much  higher  than  in 
the  low  gain  case  (Figure  4).  Note  that  the  phase  shift  at  maximum  gain  is  quite  large. 
The  gain  spectrum  is  nearly  symmetric,  and  substantial  gain  exists  at  resonance 
where  the  phase  shift  <|>  =  n/2. 

To  understand  the  mechanism  for  high  gain,  a  phase  space  simulation  using 
J  =  10°,  ao  =  1.  ar|d  v0  =  0  is  shown  in  Figure  7.  The  electrons  which  start  at 
resonance  drive  the  optical  phase  <(>,  represented  in  phase  space  by  the  shifting 

separatrix,  and  attempt  to  bunch  at  C,  =  k/2  as  discussed  earlier  when  considering 
Figure  3. 


***  FEL  Phase  Space  Evolution  *** 


3=100  a  =1  v  =0 

o  o 


Figure  7.  Phase  space  simulation  of  high  gain  FEL  with 
electrons  starting  at  resonance.  Final  gain  G  =  62  =  e4 
and  the  phase  shift  o  =  n’2.  The  separatrix  is  shown 
growing  with  the  field  strength  la  I  and  shifting  left  as  the 
phase  o  grows. 


30 


It  is  the  phase  shift  <|>  =  ji/2  which  allows  gain  as  the  relative  electron  bunch  phase 
(C  +  <(>)  =  n  drives  the  optical  field  amplitude.  This  leads  to  exponential  gain  (linear  on 
the  logarithmic  plot  of  Figure  7)  as  predicted  by  (58). 

To  summarize,  when  j  «:  1  we  must  start  the  electrons  at  v0  >  0  in  order  that 
they  bunch  near  (£  +  <|>)  =  %  and  gain  results.  For  an  FEL  with  j  »  1 ,  electrons  at 
v0  =  0  drive  the  optical  phase  with  the  high  current  until  the  bunch  is  situated  near 
(C  +  4>)  =  ft-  In  the  case  of  low  gain,  the  optical  phase  shift  is  unimportant  and  small;  in 
the  case  of  high  gain  it  is  a  crucial  aspect  of  the  interaction. 


3.  Gain  Degradation  Due  to  Electron  Beam  Quality 

In  section  C.2  of  this  chapter,  we  discussed  transverse  electron  beam  quality 
and  deduced  two  measures,  ae  and  cve,  from  the  emittance.  The  longitudinal  beam 
quality  has  a  measure 


oG  =  4nN 


A y 
Y 


(59) 


which  is  the  standard  deviation  of  a  Gaussian  spread  in  energy  (see  (17)),  and  Ay me2 
is  the  rms  energy  difference  from  the  resonant  energy  (Ref.  3). 

Consider  the  low  gain  spectrum  in  Figure  3  again.  The  full  width  of  the  main 
positive  gain  portion  at  half  the  maximum  is  approximately  n  and  the  distance  from 
maximum  gain  to  maximum  absorption  is  approximately  2k.  Therefore,  if  the  spread  in 
phase  velocities  of  the  electron  beam  Av  =  n,  we  should  expect  serious  gain 
degradation  as  many  electrons  will  be  injected  into  the  undulator  with  phase  velocity 
corresponding  to  less  than  maximum  or  even  to  negative  gain.  Inspection  of  (50a) 
and  (16)  shows  that  Av  has  contributions  from  oe,  ov4>.  and  oG. 

In  Figure  8  is  a  phase  space  simulation  using  the  same  parameters  as  those 
used  in  Figure  5  except  that  a  small  energy  spread  oG  =  1  has  been  included.  Only 
the  final  phase  space  positions  of  the  electrons  have  been  plotted.  The  electrons  are 
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spread  through  phase  space  and  bunching  has  been  degraded.  Gain  has  been 
reduced  to  G  =0.10. 


***  FEL  Phase  Space  Evolution  *** 

3=1  a  =3  V  *2.6  c_=l 

O  O  G 

0.13 


0 

0.1 

0 

-71/2  C,  3n/2  0  T  1 

Figure  8.  Phase  space  simulation  of  low  gain  FEL  with 
energy  spread  oG  =  1.  Gain  is  reduced  to  G  =  0.10. 
Compare  with  Figure  5. 


In  the  high  gain  spectrum  of  Figure  6  we  note  that  the  gain  spectrum  bandwidth 
for  /  =  100  is  approximately  2 n.  In  general,  this  width  is  proportional  to  y'1/6  in  the  high 
gain  region  (Ref.  3).  Since  there  are  no  regions  of  absorption  and  the  spectrum  is 
about  twice  as  wide  as  that  of  the  low  gain  example,  we  would  expect  the  high  gain 
FEL  to  be  more  resistant  to  beam  quality  effects.  However,  when  increasing  the 
undulator  length,  beam  quality  measures  oG,  (59),  oe,  (18),  and  ove,  (20)  increase 
proportional  to  N.  This  means  that  the  quality  of  the  electron  beam  is  magnified  by  a 
long  undulator.  Since  j  «  IN3,  where  1  is  the  peak  current,  high  gain  is  most  easily 
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and  often  obtained  by  using  a  long  undulator.  In  addition,  higher  peak  current  t  from 
the  accelerator  leads  to  poorer  beam  quality,  so  that  either  way  the  designer  tries  to 
maximize  gain  the  problem  with  beam  quality  gets  worse.  When  trying  to  get  all  one’s 
light  amplification  in  a  single  pass,  little  gain  degradation  can  be  tolerated.  Indeed  this 
may  become  one  of  the  major  drawbacks  of  using  an  FEL  amplifier  as  a  high  power 
weapon. 

In  Figure  9  is  a  phase  space  simulation  using  the  same  parameters  as  those 
used  in  Figure  7  except  in  this  case  oG  =  1.  Only  the  final  separatrix  and  electron 
phase  space  positions  have  been  plotted.  The  electrons  are  spread  throughout  phase 
space  and  bunching  has  been  degraded.  Gain  has  been  reduced  from  G  =  62  to 
G  =  50. 


***  FEL  Phase  Space  Evolution  *** 
j=100  a  =1  v  =0  G_=l 

O  O  b 


-n/2  £  3n/2  0  t  1 

Figure  9.  Phase  space  simulation  of  high  gain  FEL  with 
energy  spread  oG  =  1.  Gain  is  reduced  to  G  =  50. 
Compare  with  Figure  7,  where  G  =  62. 


33 


In  both  Figures  8  and  9  we  have  only  included  the  energy  spread  Gq 
contribution  to  beam  quality.  Clearly  the  degradation  would  be  even  worse  if  we  also 
included  the  effects  of  emittance. 

F.  SHORT  PULSES 

It  was  mentioned  during  our  derivation  of  the  wave  equation  that  the  electrons 
slip  back  relative  to  the  light  pulse  at  speed  c(  1  —  pz).  From  the  definition  of 
resonance,  we  know  that  the  distance  between  the  leading  edges  of  the  electron  and 
light  pulses  at  i  =  1  is  A/A..  This  distance  is  called  the  "slippage  length."  If  the 
electron  pulse  length  is  on  the  order  of  the  slippage  length,  the  pulse  is  considered 
short  (Ref.  3).  In  this  case  we  cannot  make  the  long  pulse  assumption  in  (44) 
because  each  section  of  the  light  pulse  sees  the  charge  density  changing  as  the 
electron  pulse  slips  past. 

In  an  optical  resonator  of  length  S,  it  takes  time  2S/c  for  the  optical  pulse  to 
make  one  round  trip.  Let  the  desynchronism  d  be  the  distance  between  the  light  and 
electron  pulses  at  x  =  0  on  each  pass  (Ref.  3).  We  normalize  d  and  the  electron 
pulse  length  oz  =  le/(NX)  to  the  slippage  length.  At  exact  synchronism,  unbunched 
electrons  slip  back  past  the  light  pulse  as  x  :  0  — >  1 ,  bunching  occurs,  but  only  the  rear 
of  the  light  pulse  gets  amplified.  On  successive  passes  the  delay  in  bunching 
becomes  longer  and  the  amplification  moves  further  back  on  the  light  pulse.  Since 
only  the  rear  of  the  light  pulse  is  being  amplified  while  the  losses  (as  determined  by 
the  resonator  quality  factor  Q)  are  the  same  at  all  points  in  the  pulse,  its  centroid 
moves  at  a  speed  slower  than  c.  This  is  called  "lethargy,"  and  ultimately  will  lead  the 
pulse  out  of  the  region  of  gain  (that  is,  it  won’t  overlap  the  electron  pulse)  and  the 
pulse  amplitude  will  vanish  (Ref.  3). 

To  compensate  for  the  effects  of  lethargy  a  small  amount  of  desynchronism 
must  be  added.  This  is  accomplished  by  reducing  the  path  length  S  that  the  light 
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pulse  must  travel  by  adjusting  the  length  of  the  optical  cavity.  In  this  case, 
d  =  -AS /(NX)  (Ref.  3). 

Simulations  of  an  FEL  using  picosecond  electron  pulses  will  be  discussed  in 
Chapter  III.  There  we  will  explore  the  effects  of  desynchronism  on  the  power 
evolution. 


G.  FEL  WITH  GAUSSIAN  OPTICAL  MODE 


In  section  E  on  gain  we  discussed  the  interaction  of  the  electron  beam  with 
light  approximated  by  plane  waves.  An  optical  cavity  bounded  by  spherical  mirrors 
supports  Hermite-Gaussian  (and  other)  modes,  which  are  solutions  of  (45)  instead  of 
(47)  and  accurately  describe  the  laser  beam  inside  and  outside  the  cavity  (Ref.  7). 
The  optical  beam  can  be  constructed  from  any  number  of  these  modes,  but  the  one 
most  desired  in  a  weapons  laser  for  its  propagation  characteristics  is  called  the 
fundamental  Gaussian  mode,  or  just  Gaussian  mode. 

The  usual  Gaussian  mode  as  given  in  Ref.  7  can  be  written  in  terms  of  our  FEL 
dimensionless  parameters.  The  result  is 


where 


=  _^e-rWroe<'0G(t) 

W(t) 


<1 >g(x)  =  -tan 


X  -  Xu 


-  xw) 
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(61) 


is  the  phase  shift  of  the  Gaussian  mode  relative  to  that  of  a  plane  wave, 


w^x)  =  1  + 


x  -  x„ 


l  z°  J 


(62) 


is  the  square  of  the  mode  radius,  which  is  normalized  to  the  mode  waist  radius  , 
and 
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zo  =  wo  (63) 

is  the  dimensionless  Rayleigh  length  (Ref.  3).  The  mode  waist  position  is  xw.  Radial 
dimensions  are  normalized  to  VlA/jc,  the  characteristic  mode  size.  We  see  in  (62)  that 
zo  is  the  distance  over  which  the  mode  area  doubles.  In  a  symmetric  optical  cavity 
with  spherical  mirrors  a  distance  S  apart  and  radii  of  curvature  R,  the  dimensionless 
Rayleigh  length  is  given  by  z0  =  (2 SR  -  S2)'/2/2L  where  L  is  again  the  undulator 
length  (Ref.  7).  The  Rayleigh  length  Zq  and  the  waist  position  xw  completely 
characterize  the  Gaussian  mode  in  the  optical  cavity. 

The  amplitude  of  the  Gaussian  mode  la  I  is  a  maximum  at  x  =  xw  where 
w2(t)  =  1.  It  decreases  away  from  the  waist  in  both  directions  and  is  a  minimum  at 
both  t  =  0  and  t  =  1  (if  xw  =  0.5).  At  any  time  x,  the  amplitude  decreases  in  the  radial 
direction  as  a  Gaussian,  hence  the  name.  The  phase  of  the  Gaussian  mode  <{>G  (61) 
contains  two  terms.  The  first  term  provides  a  phase  shift,  called  the  Guoy  phase  shift, 
relative  to  the  plane  wave,  and  switches  sign  at  the  waist  (Ref.  7).  The  second  term  in 
(61)  decreases  the  phase  shift  for  positions  off-axis  which  creates  the  wavefront 
curvature  necessary  for  focussing  the  beam  at  the  waist. 

In  the  electron  phase  space,  the  term  ao/w(x)  causes  the  height  of  the 
separatrix  to  increase  as  x  :  0  xw  and  then  decrease  as  x  :  xw  1 .  The  effect  is 
proportional  to  w“,/j(x)  (from  our  knowledge  of  the  separatrix  height)  and  complicates 
the  phase  space  paths.  The  phase  of  the  Gaussian  mode  (61)  is  continuously  shifted 
as  x  :  0  1 .  The  shift  has  the  sign  opposite  to  that  caused  by  the  FEL  interaction  in 

(48b).  It  is  therefore  represented  in  phase  space  by  the  separatrix  shifting  to  the  right 
as  x  :  0  -»  1 .  Looking  at  Figure  2  again,  we  can  predict  that  the  electrons  will  have  to 
start  at  a  phase  velocity  v0  >  2.6  in  order  to  compensate  for  the  shift  and  still  bunch  at 
the  correct  phase  for  gain.  (Compare  this  to  the  case  of  the  high  gain  FEL  where  the 
phase  velocity  had  to  be  lowered  in  order  to  bunch  at  the  correct  phase.) 
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In  order  to  observe  the  effects  of  the  Gaussian  mode  on  the  FEL  interaction, 
assume  z0  »  1,  nearly  plane  waves,  and  expand  (61)  and  (62)  to  first  order  in  (1/z0). 
The  phase  becomes  <|>G  =  -(t  -  iw)lz0,  and  w(x)  ~  1.  The  pendulum  equation  (36)  in 

oo 

the  presence  of  the  Gaussian  mode  is  C,  =  I aG  lcos(£  +  <J>G)  which  for  z0  »  1  becomes 


_  -r2/z0 

=  aQe  0  cos 

Co  +  — + 

c  ■> 

V.--S- 

X 

i 

N 

o 

t - 

N 

o 

\ _ 

(64) 


when  expanded  in  weak  fields  to  first  order  in  a0  (Ref.  3).  Comparing  (64)  with  (36) 

oo 

we  see  three  distinct  differences.  First,  the  amplitude  of  the  phase  acceleration  £ 
decreases  for  electrons  off  axis  because  they  see  a  weaker  field  strength 
la  I  -»  a0exp(-r2/z0).  Second,  there  is  a  shift  in  the  initial  phase  C,0  of  all  electrons  in 
the  amount  xw/z0.  This  is  inconsequential  since  the  electrons  are  initially  distributed 
uniformly  throughout  phase  space.  The  last  difference  is  in  the  phase  velocity,  which 
is  decreased  by  Av  =  1/z0.  This  leads  to  a  wavelength  change  of  Ak/X=  V(2nNz0) 
from  (50b).  This  also  requires  an  increase  in  the  initial  phase  velocity  for  maximum 
gain  to  v™3*  =  2.6  +  1/z0  in  the  low  gain  FEL.  This  effect  is  caused  by  the  optical 
phase  shift  (61)  and  is  shown  in  Figure  10.  The  figure  is  the  result  of  a  simulation 
similar  to  that  in  Figures  4  and  6  except  that  the  finite  size  of  the  optical  mode  is 
included.  Here  the  Rayleigh  length  z0  =  5  and  maximum  gain  is  found  at 
v0  =  2.9  =  (2.6  +  1/5). 


The  presence  of  the  Gaussian  optical  mode  produces  measurable  effects  in  the 
performance  of  the  FEL.  This  is  the  theme  of  the  bulk  of  this  work,  including  the  final 
chapter  on  Compton  backscattering  where  we  calculate  the  photon  emission  rate  for 
electrons  in  the  Gaussian  field  of  a  laser. 


37 


***  FEL  3d  simulation,  single-pass  gain  *** 
3=1  zo~5  ao=0.01 
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0.0' 


0.0 
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Figure  10.  Gain  spectrum  for  low  gain  FEL  with  Rayleigh 
length  z0  =  5.  Maximum  gain  is  achieved  at  v0  =  2.9, 
shifted  from  v0  =  2.6  by  the  Gaussian  mode. 


38 


III.  SIMULATIONS  OF  FIREFLY  1KW  FEL 


A.  BACKGROUND 

The  FIREFLY  (Far  Infrared  FEL)  IkW  FEL  at  Stanford  University  was  proposed 
in  1994  as  an  upgrade  to  the  existing  Infrared  FEL  at  Stanford  University  (Ref.  8). 
With  four  new  10  MeV  accelerator  structures  and  a  5  MeV  electron  injector,  an 
electron  beam  with  energy  45  MeV  could  be  generated.  The  electron  injector 
operating  at  1 1 .8  MHz  would  produce  an  average  current  of  T  =  1  ma,  and  the 
average  power  in  the  electron  beam  would  be  Pe  =  TV  =  45  kW.  A  good  measure  of 
the  maximum  efficiency  of  an  FEL,  i.e.,  how  much  energy  one  can  extract,  is  the  gain 
spectrum  bandwidth.  Considering  Figure  4  and  (50b),  we  find  that  Ay/y  =  1/(2 N).  If 
we  assume  that  we  can  extract  the  full  Ay/y=  1/(2 N)  from  the  beam  with  N  =  25,  then 
the  average  power  out  of  the  cavity  would  be  P  =  Pe/(2N)  =  900  W.  This  average 
power  would  be  almost  1 00  times  greater  than  the  world’s  most  powerful  FEL,  and  it 
was  estimated  this  could  be  done  for  only  one  million  dollars  with  existing  technology. 

The  proposed  FEL  and  two  variants  with  lower  electron  beam  energies  of  25 
MeV  and  1 0  MeV  were  simulated.  All  other  parameters  were  kept  the  same  and  are 
listed  in  Table  1. 

B.  SIMULATIONS  AND  RESULTS 

Because  all  three  FELs  simulated  have  emission  wavelengths  X  in  the  far 
infrared,  the  characteristic  mode  size  VLX/it  is  large  enough  that  diffraction  (i.e.,  the 
Gaussian  mode)  should  not  dominate  the  interaction.  The  slippage  length  NX, 
compared  to  the  electron  pulse  length,  is  also  large,  so  the  effects  of  short  optical 
pulses  must  be  addressed.  (Some  basic  consideration  of  short  pulse  effects  was  done 
in  section  F  of  Chapter  II.)  To  study  short  pulse  effects,  simulations  following 
longitudinal  multimodes  were  used.  These  simulations  divide  the  optical  pulse  into  a 
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number  of  fixed  longitudinal  field  sites,  or  bins,  a(z)  (Ref.  3).  The  simulations  are 
called  longitudinal  multimode"  because  dividing  the  optical  pulse  into  longitudinal 
sites  a  ->  a(z)  is  equivalent  to  a  -»  a(kz)  in  Fourier  space;  the  number  of  sites  z  is 
equivalent  to  a  number  of  modes  kz. 


1 

II 

III 

Beam  Energy 

45  MeV 

25  MeV 

10  MeV 

Peak  Current 

40  A 

40  A 

40  A 

Micropulse  Length 

2.1  ps 

2.1  ps 

2.1  ps 

Undulator  Type 

Linear 

Linear 

Linear 

K 

1.1 

1.1 

1.1 

N 

25 

25 

25 

Xq 

6.0  cm 

6.0  cm 

6.0  cm 

X 

8.4  pm 

26.6  pm 

156.7  pm 

j 

1.7 

2.9 

7.2 

3.2 

1.0 

0.2 

Table  1.  Parameters  for  FIREFLY  IkW  FEL  (column  I) 
and  two  variants  (columns  II  and  III)  used  in  simulations. 


The  slower  moving  electrons  slip  back  past  the  optical  pulse  to  sites  z  -  x  and  interact 
with  a  number  of  field  sites  corresponding  to  the  slippage  distance  NX  (recall  that  z  is 
normalized  to  the  slippage  distance).  Conversely,  each  field  site  a  (z)  interacts  with  a 
length  NX  along  the  electron  micropulse  as  the  micropulse  slips  back.  We  generalize 
the  pendulum  (36)  and  wave  (47)  equations  to  include  the  slippage: 
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=  lazlcos(k+T  +  <|>z) 

(65) 

:  ~Jz-n<e  M>z+t  • 

(66) 

As  t :  0  ->  1  in  increments  of  Ax  in  the  simulation,  (65)  and  (66)  show  how  the  field 
changes  at  site  z  are  caused  by  the  electron  dynamics  at  the  micropulse  site  z  +  s  Ax, 
where  s  =  0,  1,  2,  ...,Nt  is  an  integer  and  Nt  =  1/Ax  is  the  total  number  of  time  steps, 
and  vice  versa. 

Contained  in  Figure  1 1  is  the  output  from  a  simulation  of  the  1  kW  machine, 
whose  parameters  are  listed  in  column  I  of  Table  1.  The  dimensionless  current 
density  is  j  =  1 .7,  the  micropulse  length  is  az  =  3.2,  cavity  losses  are  determined  by 
0  =  100,  and  the  desynchronism  d  =  0.03  is  large  (see  Chapter  II  section  F  for 
definition  of  d). 


****  FEL  Pulse  Evolution  **** 

j=1.7  C  =3.2  Q=100  d=0 . 03 

z 


-20  *  .  20  -19  v  19  0  n  2000 

Figure  11.  Longitudinal  multimode  simulation  of 


FIREFLY. 
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The  parabolic  electron  pulse  is  drawn  in  the  lower  left  window  at  x  =  0  and  at  %  =  1  to 
illustrate  slippage.  Above  the  electron  pulse  is  plotted  the  optical  field  amplitude 
each  pass  n.  The  scale  ranges  from  zero  in  black  to  the  maximum  in 
white;  a  contour  line  is  drawn  at  the  median  field  amplitude.  Above  the  field  evolution 
is  plotted  the  final  optical  field  amplitude  la(z,n  =  2000)1  in  black.  At  the  bottom 
center  is  drawn  the  gain  spectrum  (57).  Above  this  is  the  optical  power  spectrum 
evolution  P(v,n),  with  the  final  power  spectrum  P(v,n  =  2000)  drawn  on  top  in  black. 
The  power  spectrum  P(v,n)  is  the  Fourier  transform  of  the  square  of  the  optical  field 
la(z,n)l2  at  each  pass  n.  The  phase  velocity  v  in  the  power  spectrum  is  used  in  the 
sense  of  (50b);  that  is,  it  represents  the  wavelength  X  so  that  P(v,n)  shows  the 
linewidth  of  the  radiation.  The  total  power  evolution  P(n)  is  plotted  at  lower  right  in 
the  figure.  Above  this  is  plotted  the  electron  phase  velocity  distribution  f(v,n).  The 
distribution  f(v,n)  is  a  projection  of  the  final  electron  phase  space  positions  (£,v)  onto 
the  v-axis  at  each  pass  n .  As  such,  v  in  this  plot  is  to  be  understood  in  the  sense  of 
(50a);  that  is,  it  represents  the  electron  energy  ymc2  so  that  f  (v,n)  is  interpreted  as  an 
energy  spread. 

Having  described  all  the  plots  in  Figure  11,  we  now  interpret  their  features. 
The  large  desynchronism  has  caused  the  optical  pulse  la(z,n)l  to  be  artificially  long 
by  advancing  the  front  edge  of  the  pulse  ahead  of  the  electrons.  The  electron  pulse 
amplifies  only  the  trailing  edge  of  the  optical  pulse  on  each  pass.  This  leads  to 
reduced  coupling;  the  peak  field  strength  is  only  la  I  =  25.  As  evidenced  by  the  power 
spectrum  P(v,n),  the  long,  smooth  optical  pulse  leads  to  a  narrow  linewidth.  The 
electron  phase  velocity  distribution  f(v,n)  does  not  evolve  any  further  after  the  first 
few  hundred  passes.  At  saturation,  the  electrons  reach  the  bottom  of  the  separathx  in 
phase  space  (Ref.  3).  This  means  that  the  electron  distribution  f(\\n)  should  have  a 
width  of  about  4VTaT  =  18  at  saturation,  the  peak-to-peak  height  of  the  separatrix  at 
la  I  =20,  the  final  field  amplitude.  In  the  figure,  the  width  of  f(v,n)  is  about  18, 
showing  good  agreement.  The  power  P(n)  at  each  pass  n  is  obtained  by  integrating 
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over  the  optical  pulse  la(z,n)l2.  The  power  evolution  P(n)  plotted  at  lower  right 
shows  that  a  steady  state  is  reached  after  approximately  n  =  700  passes  where 
P{n)  =  91. 

Plotted  in  Figure  12  are  the  results  of  several  simulations  like  that  above.  The 
steady  state  dimensionless  power  is  plotted  versus  the  desynchronism. 


Figure  12.  Plot  of  outcoupled  dimensionless  power 
versus  desynchronism  for  FIREFLY. 


The  maximum  on  the  curve  at  d  =  0.0005  corresponds  to  an  actual  output  average 
power  P  =  900  W,  the  desired  output,  but  requires  the  injector  to  operate  at  23.6  MHz, 
or  double  the  rate  discussed  earlier.  Evidently  the  short  pulses  lead  to  reduced 
coupling  between  the  electrons  and  optical  field,  lowering  FEL  efficiency.  Doubling  the 
operating  frequency  of  the  electron  injector  was  proposed  in  Ref.  8  as  a  follow-on  to 
the  1  kW  upgrade  for  the  purpose  of  attaining  2  kW  average  power.  Note  that  the 
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maximum  in  Figure  12  is  rather  sharp.  Since  the  slightest  change  in  of  can  produce  a 
large  change  in  the  output,  the  FEL  is  less  stable  when  operated  in  this  region  (Ref. 
3). 

The  FIREFLY  variant  listed  in  column  II  of  Table  1  has  a  shorter  dimensionless 
electron  pulse  length  oz  than  that  in  Figure  11.  The  output  of  a  simulation  with 
j  =  2.9,  cz  =  1,  Q  =  100,  and  d  =  0.005  is  shown  in  Figure  13.  The  features  of  the 
plots  show  the  occurrence  of  the  trapped-particle  instability  (Ref.  3). 


****  FEL  Pulse  Evolution -  ***■* 
j=2.9  0Z=1  Q=100  d«=0.005 


-4  z  4  -25  v  25  0  n  2000 

Figure  13.  Longitudinal  multimode  simulation  for  the 
variant  of  FIREFLY  listed  in  column  II  of  Table  1. 


This  instability  occurs  when  the  field  amplitude  la  I  becomes  so  strong  that  electrons  in 
phase  space  traverse  their  entire  path  (see  Figure  2)  during  a  single  pass  through  the 
undulator.  executing  a  synchrotron  oscillation  with  frequency  vs  =  VTaT  =  2n  (Ref.  3). 
In  real  space,  this  oscillation  represents  motion  superimposed  on  the  electron 
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trajectory  and  leads  to  emission  into  sidebands  v0  ±  vs  (Ref.  3).  As  these  sidebands 
grow,  the  FEL  power  will  increase  and  the  fields  will  get  stronger,  causing  more 
electrons  to  execute  synchrotron  oscillations,  etc.  At  least  four  sidebands  can  be  seen 
in  the  power  spectrum  P(v,n)  in  Figure  13.  These  sidebands  modulate  the  optical 
pulse  envelope  la(z,n)l  as  shown,  and  the  increased  field  strength  causes  the  wide 
electron  distribution  f(v,n).  In  this  case,  the  instability  has  imposed  a  periodic 
structure  on  the  optical  field  envelope  \a(z,n)\  and  the  power  spectrum  P(v,n)  which 
results  in  oscillations  in  the  total  power  P(n)  over  many  passes  n. 

Figure  14  plots  the  dimensionless  power  versus  desynchronism  from  several 
simulations  like  that  of  Figure  13. 


Figure  14.  Plot  of  dimensionless  outcoupled  power 
versus  desynchronism  for  variant  of  FIREFLY  listed  in 
column  II  of  Table  1. 
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We  see  from  the  plot  that  Figure  12  with  d  =  0.005  lies  in  the  transition  region  in  the 
desynchronism  curve,  where  the  fields  are  strong  enough  to  cause  the  trapped-particle 
instability,  but  where  operation  of  the  FEL  is  relatively  stable.  The  peak  in  this  curve 
at  d  =  0.003  corresponds  to  an  outcoupled  average  power  of  P  =  140  W. 

To  further  illustrate  the  effects  of  slippage,  the  example  listed  in  column  III 
represents  FIREFLY  operating  at  a  much  longer  wavelength  X  =  156.7  pm.  The 
slippage  length  is  therefore  long  and  leads  to  a  small  value  of  az  =  0.2.  Figure  15 
contains  the  result  of  a  simulation  using  j  =  7.2,  az  =  0.2,  Q  =  100,  and  the  very  small 
value  of  d  =  0.0005. 


****  FEL  Pulse  Evolution  **** 

3=7.2  Oz=0.2  Q=100  d=0. 0005 


| a(z,n)  |  0—133 _ P(v,n)  f(v,n) 


-4  z  4  -101  v  101  -101  v  101 


-4  z  4  -101  v  101  0  n  6000 


Figure  15.  Longitudinal  multimode  simulation  for  vanant 
of  FIREFLY  listed  in  column  III  of  Table  1. 


The  small  desynchronism  means  that,  unlike  in  Figure  11.  the  light  pulse  a(z,n) I  has 
not  been  artificially  lengthened.  The  optical  pulse  reaches  steady  state  with  a  length 
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comparable  to  the  electron  pulse  length  shown  in  the  plot  at  lower  left.  The  fields  are 
strong  and  the  trapped-particle  instability  again  leads  to  structure  in  the  optical  pulse. 
We  should  note  here  that  the  spikes  to  the  left  of  the  main  optical  pulse  are  only  about 
one  wavelength  of  light  in  length.  Because  this  would  violate  the  slowly-varying 
amplitude  assumption  in  the  wave  equation  (47),  these  may  be  meaningless.  The 
short,  modulated  optical  pulse  results  in  a  broad  power  spectrum  P(v,n)  (note  the 
scale  in  v  compared  with  that  in  Figure  13).  The  strong  fields  also  result  in  a  broad 
electron  distribution.  Maximum  power  P(n)  =  900  is  achieved  at  n  ~  1500  passes,  but 
steady-state  P(n)  =  520  isn’t  reached  until  n  -  4000  passes  which  is  predicted  to  take 
about  0.3  ms  in  FIREFLY. 

Figure  16  plots  power  versus  desynchronism  for  column  III  of  Table  1. 


Figure  16.  Plot  of  outcoupled  dimensionless  power 
versus  desynchronism  for  variant  of  FIREFLY  listed  in 
column  III  of  Table  1. 
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We  note  that  d  =  0.0005,  that  used  in  Figure  15,  is  the  location  of  the  peak  in  this 
curve  and  corresponds  to  an  average  power  of  P  =  17  W. 

C.  CONCLUSIONS 

The  operational  characteristics  of  far  infrared  wavelength  FELs  driven  by  short, 
picosecond  electron  micropulses  will  be  determined  by  short  pulse  effects.  The 
desynchronism  coupled  with  the  cavity  Q  will  affect  the  length  and  structure  of  the 
optical  pulse.  This  in  turn  determines  the  width  and  structure  of  the  power  spectrum. 

Based  on  our  simulations  of  the  1  kW  FEL,  as  depicted  in  Figures  11  and  12, 
the  goal  of  1  kW  average  power  can  be  achieved,  but  at  double  the  electron  injector 
duty  cycle  first  proposed.  More  simulations  should  be  done  to  substantiate  this 
possibility.  Doubling  the  injector  duty  cycle  was  proposed  in  Ref.  8,  but  as  a  follow-on 
to  the  1  kW  proposal  as  a  means  of  achieving  2  kW.  We  also  note  that  the  trapped- 
particle  instability  will  probably  occur. 
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IV.  TRANSVERSE  EFFECTS  IN  UV  FREE  ELECTRON  LASERS 


A.  BACKGROUND 

In  an  ultraviolet  free  electron  laser  (UVFEL),  the  characteristic  optical  mode 
size  VH/7C  tends  to  be  small  because  of  the  short  optical  wavelength  X.  For  example, 
if  X  =  200  nm  and  L  =  3  m,  VTX/tc  =  0.4  mm.  This  is  also  approximately  the  size  of 
the  electron  beam.  For  the  ultraviolet  FEL,  the  effects  of  the  spatial  extent  of  the 
electron  beam  as  compared  to  the  optical  mode  are  of  concern. 

The  UVFEL  under  consideration  in  this  chapter  is  that  proposed  by  the  Laser 
Processing  Consortium  (LPC)  at  Thomas  Jefferson  National  Accelerator  Facility 
(TJNAF)  in  Newport  News,  Virginia  (Ref.  9).  The  LPC  is  a  group  of  industry  and 
university  scientists  and  engineers  interested  in  an  efficient,  reliable,  coherent  UV  light 
source  for  processing  materials.  On  the  order  of  100  kW  of  UV  light  is-  needed  to 
process  commercial  products  ranging  from  nylon  to  steel.  The  parameters  for  the 
UVFEL  proposed  in  Ref.  9  are  listed  in  Table  2. 


Beam  Energy 

200  MeV 

Peak  Current 

270  A 

Micropulse  Length 

0.25  ps 

Undulator  Type 

Linear 

K 

0.95 

N 

72 

*0 

3.3  cm 

X 

204  nm 

Table  2.  Parameters  for  the  LPC  proposed  UVFEL. 
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In  this  chapter  we  analyze  the  LPC  UVFEL  using  transverse  multimode 
simulations.  These  simulations  are  described  in  section  B.  We  use  these  simulations 
in  section  C  to  study  the  effects  of  the  optical  mode  size  on  FEL  gain.  In  section  D  we 
study  the  effects  of  electron  beam  focussing  on  FEL  gain. 


B.  TRANSVERSE  MULTIMODE  SIMULATIONS 


In  Chapter  III  we  extended  the  pendulum  and  wave  equations  in  the  z-direction 
in  order  to  properly  simulate  longitudinal  multimodes  by  dividing  the  optical  pulse  into 
■z  sites.  To  simulate  transverse  multimodes,  we  make  the  same  transition  in  the  x- 
and  y-directions.  The  optical  field  is  extended  to  sites  in  x  and  y  so  that  a  -»  a(x,y). 
Generalizing  the  wave  equation  (45)  and  the  pendulum  equation  (36)  gives  (Ref.  3): 


4v| + ± 

4  +  8x 


a(x,y,x)  =  -<y  exp(-/;(x,y,x))>(xyT) 


C  (*,y,x)  = 


a(x,y,x) 


cos 


£(x,y,x)  +  <t>(x,y,x) 


(67) 


(68) 


The  transverse  dimensions  are  normalized  to  the  characteristic  mode  size  VXDtc.  The 
set  of  equations  (67)  and  (68)  are  fully  self-consistent  in  three  dimensions  (x,y,x). 


For  the  UVFEL  represented  in  Table  2,  we  calculate  an  electron  pulse  length 
slightly  greater  than  five  slippage  distances  a2  =  5.1.  Therefore,  it  can  be  assumed 
that  slippage  will  not  be  an  important  consideration,  so  that  we  may  follow  only  one 
representative  longitudinal  site  as  in  (67)  and  (68). 

In  section  G  of  Chapter  II,  we  discussed  the  fundamental  Gaussian  mode  but 
mentioned  that  in  general  any  number  of  Hermite-Gaussian  transverse  modes  are 
present.  In  order  to  model  the  evolution  of  the  FEL  from  weak  fields  through 
saturation,  these  higher-order  modes  are  included  as  the  FEL  interaction  amplifies  and 
distorts  the  optical  mode.  Both  the  optical  field  amplitude  and  phase  can  be  altered 
(Ref.  3).  An  example  of  this  is  shown  in  Figure  17. 
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Figure  17  shows  the  evolution  of  the  FEL  from  weak  fields  a0  =  0.0001  through 
saturation  using  j  =  20  and  a  small  energy  spread  oG  =  0.5.  The  Rayleigh  length  for 
the  cavity  is  z0  =  0.5.  The  top  set  of  plots  show  the  FEL  oscillator  after  n  =  12 
passes.  The  electron  phase  space  plot  shows  no  noticeable  bunching  and  the  final 
gain  G  =  0.94  is  in  steady  state  (i.e.,  G(n)  constant  )  with  power  increasing 
exponentially.  At  right  is  shown  the  optical  field  amplitude  between  the  mirrors 
la(x,  t)I.  The  mode  is  clearly  distorted  from  a  fundamental  Gaussian  mode  with  the 
peak  field  (in  black)  located  at  x  =  1  where  single  pass  gain  is  the  maximum.  Moving 
down  the  set  of  figures,  the  electron  phase  space  begins  to  exhibit  bunching  at 
n  =  24.  The  bunch  descends  along  the  v-axis,  indicating  energy  loss,  and  near 
saturation  at  n  =  28  it  starts  to  ascend  taking  energy  back  from  the  optical  field.  At 
saturation  (n  =  40),  the  electron  bunch  is  clearly  moving  up  in  phase  space  at  x=  1. 
The  motion  of  the  electron  bunch  through  phase  space,  caused  by  the  strengthening 
fields,  is  attended  by  a  drop  in  final  gain  G(x=1)  as  the  electrons  are  absorbing  light  at 
the  end  of  the  undulator.  This  means  that  on  each  pass,  G(x)  peaks  at  x  <  1  along 
the  undulator  and  the  peak  field  moves  back  along  the  undulator.  Scanning  down  the 
page  through  the  plots  of  la(x,x)l,  we  see  the  peak  field  strength  moving  back  along 
the  undulator.  At  saturation,  the  peak  gain  per  pass  is  nearly  at  the  middle  of  the 
undulator  and  the  final  gain  is  reduced  to  G  =0.0034.  The  optical  field  la(x,  x)l 
resembles  the  fundamental  Gaussian  mode  at  saturation  where  the  gain  is  low. 

The  evolution  of  the  distorted  optical  mode  to  the  Gaussian  requires  the  use  of 
the  transverse  multimode  simulations  described  above.  Multiple  transverse  modes  are 
included  in  the  simulation  to  construct  the  distorted  optical  fields  of  Figure  17  for  the 
same  reason  it  takes  a  large  number  of  sine  waves  to  synthesize  a  square  pulse  in 
Fourier  analysis. 
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c. 


OPTICAL  MODE  SIZE  EFFECT  ON  FEL  GAIN 


It  has  been  shown  that  the  Rayleigh  length  which  optimizes  gain  for  a  filament 
electron  beam  (i.e.,  oe  =0)  is  z0  =  1/Vi2  =  0.3  (Refs.  3,  10).  This  is  the  Rayleigh 
length  which  maximizes  the  average  filling  factor  F  (averaged  along  the  undulator 
length)  by  minimizing  the  volume  occupied  by  the  Gaussian  optical  mode  around  the 
electron  beam.  Reference  10  also  discusses  gain  in  a  Gaussian  mode  for  non-zero 
electron  beam  radius,  but  again  using  arguments  based  on  the  filling  factor.  The 
conclusion,  not  surprisingly,  is  that  maximum  gain  is  achieved  when  the  Rayleigh 
length  is  chosen  such  that  the  electron  beam  fits  inside  the  optical  mode  (recall 
zo  =  wo  )•  But  in  the  case  of  the  UVFEL,  it  may  be  necessary  to  have  an  electron 
beam  bigger  than  the  optical  mode  because  of  the  small  mode  size  'JXLJk.  In  that 
case,  one  cannot  quantify  gain  in  terms  of  the  filling  factor,  which  may  be  greater  than 
one.  The  three  dimensional  simulations  described  in  the  previous  section  are  used  to 
find  the  influence  of  the  Rayleigh  length  on  FEL  gain. 

Figure  18  shows  the  result  of  a  three-dimensional,  transverse  multimode 
simulation  of  the  UVFEL  parameters  listed  in  Table  2.  The  corresponding 
dimensionless  parameters  are  current  density  j  =  30  with  F  is  no  longer  included  in  j, 
v0  =  4.5  is  the  optimum  initial  phase  velocity,  the  Rayleigh  length  is  z0  =  0.3, 
a0  =  0.001  is  the  initial  optical  field  amplitude,  the  cavity  loss  is  Q  =  5,  oe  -  0.45  is  the 
dimensionless  electron  beam  radius,  oe  =  0.9  is  the  angular  spread,  the  energy  spread 
is  cG  =  1.3,  and  top  =  1.0  is  the  betatron  frequency.  The  dimensionless  emittance 
e  =  0.43  corresponds  to  an  actual  normalized  emittance  £W  =  11  jr  mm-mrad 
(c*  =  y^rms)-  Tb®  plot  at  top  left  shows  the  evolution  of  the  optical  field  amplitude  at 
the  end  of  the  undulator  with  each  pass  la(x,n)l.  To  the  right  is  plotted  the  final 
optical  field  amplitude  at  the  right-hand  mirror  la(x,y)l.  The  optical  field  between  the 
mirrors  la(x,t)l  at  y  =  0  is  plotted  in  the  center;  the  tick  marks  at  t  =  0  and  t  =  1  show 
the  beginning  and  end  of  the  undulator.  Overlaid  on  this  plot  are  the  positions  of  a 
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random  sample  of  the  electrons  used  in  the  simulation.  In  this  figure,  oe/w0  =  0.8  and 
the  electrons  near  the  edges  of  the  Gaussian  mode  see  weaker  optical  fields  and  thus 
a  degraded  interaction  (see  (60)). 

TRANSVERSE  PEL  MODEL 


Figure  18.  Transverse  multimode  simulation  of  UVFEL. 
Gain  is  degraded  due  to  large  electron  beam 
radius/optical  mode  waist  radius  ratio  ce/w0  =  0.8. 


The  plots  along  bottom  show  the  electron  phase  velocity  distribution  at  each  pass 
the  final  phase  space  positions  (£,v)  of  a  sample  of  electrons,  and  the  power 
P(n)  and  gain  ln(1+G(n))  evolutions.  The  gain  G(n)  is  plotted  at  the  end  of  the 
undulator  G(t=1)  after  each  pass  n.  These  plots  show  FEL  gain  in  steady  state  (i.e., 
G(n)  approximately  constant)  with  final  gain  G  =  0.75  and  power  increasing 
exponentially.  Electrons  in  phase  space  are  spread  out  due  to  oG  and  o0  so  that 
bunching  can  barely  be  seen  at  £  =  n  even  though  the  limit  of  weak  fields  la  I  =  n  has 
been  reached. 
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The  results  of  many  simulations  like  that  in  Figure  18  were  used  to  plot  gain 
normalized  to  the  maximum  one-dimensional  gain,  G0  =  0.135/F,  where 
F  =  KGg/{z0  +  1/12z0)  is  the  average  filling  factor,  versus  Rayleigh  length  (Refs.  3,  10). 
The  result  is  plotted  in  Figure  19. 


Figure  19.  Plot  of  normalized  gain  versus  Rayleigh 
length  for  UVFEL.  Normalized  gain  is  minimum  where 
optical  mode  volume  is  minimum  z0  =  0.3. 


The  first  feature  to  note  in  the  figure  is  that  for  the  UVFEL,  z0  =  0.3  corresponds  to  the 
minimum  normalized  gain  because  the  optical  mode  volume  is  minimized.  The  optical 
mode  volume  must  be  large  enough  to  envelope  the  electron  beam.  Decreasing  the 
Rayleigh  length  below  z0  =  0.3  causes  the  mode  to  expand  quickly  away  from  the 
mode  waist.  This  makes  the  mode  large  at  the  ends  of  the  undulator,  enveloping  the 
electron  beam,  and  increases  normalized  gain.  Increasing  the  Rayleigh  length  from 
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Zq  —  0.3  makes  the  optical  mode  larger  everywhere  and  increases  gain  until  the  mode 
gets  so  large,  z0  ~  0.9  in  this  case,  that  electrons  see  weaker  fields  everywhere  and 
gain  is  reduced.  For  a  larger  electron  beam,  the  optimum  value  of  z0  for  normalized 
gain  would  be  larger  than  z0  =  0.9. 

D.  ELECTRON  BEAM  SIZE  EFFECT  ON  FEL  GAIN 

Now,  given  that  the  Rayleigh  length  of  the  optical  cavity  has  been  chosen,  how 
is  FEL  gain  affected  by  the  choice  of  electron  beam  radius?  Since  cop  =  1 .0  for  this 
UVFEL,  the  electron’s  betatron  trajectories  are  nearly  straight  lines,  allowing  them  to 
be  focussed  at  the  undulator  center  using  external  magnets.  Keeping  emittance 
e  =  oe-\[c^  fixed,  decreasing  the  electron  beam  radius  at  the  waist  increases  the 
angular  spread  and  degrades  the  interaction,  so  there  must  be  a  limit  to  how  small 
one  can  focus  the  electron  beam.  In  all  simulations,  the  Rayleigh  length  is  held  fixed 
at  z0  =  0.5. 

Figure  20  contains  the  output  of  a  simulation  of  the  UVFEL  electron  beam  with 
e  =  0.43.  In  this  example,  we  have  reduced  the  electron  beam  radius  to  oe  =  0.3, 
which  increases  angular  spread  to  oe  =  2.1.  We  also  were  required  to  increase  the 
current  density  to  j  =  68  in  order  to  keep  peak  current  constant.  The  optimum  initial 
phase  velocity  is  v0  =  4.3  and  all  other  parameters  are  the  same  as  those  in  Figure 
18.  The  phase  space  plot  (£,v)  is  dominated  by  the  (exponential)  angular  spread 
which  indicates  ae  has  been  made  too  small.  In  the  plot  of  la(x,x)l  in  the  center,  one 
can  see  the  electrons  focussed  to  ae  =  0.3  at  the  center  of  the  undulator.  but  the  large 
angular  spread  cames  many  electrons  outside  the  optical  mode  at  t  =  0  and  t  =  1 . 
This  leads  to  gain  degradation. 

Figure  21  plots  the  normalized  gain  versus  electron  beam  radius  for  three 
different  omittances:  (a)  c  =  0.28.  (b)  e  =  0.43,  and  (c)  e  =  0.6.  (The  simulation  in 
Figure  20  is  from  curve  (b).)  We  see  in  the  curves  that  decreasing  the  electron  beam 
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radius  from  ae  =  1 .0  increases  gain  as  more  electrons  are  enveloped  in  the  optical 
mode.  The  beam  radius  can  be  made  too  small,  however,  causing  gain  to  decrease 
because  of  the  increasing  angular  spread. 


TRANSVERSE  PEL  MODEL 


Figure  20.  Transverse  multimode  simulation  of  UVFEL. 
Electron  beam  focussed  at  x  =  0.5  but  angular  spread 
causes  beam  to  expand  outside  optical  mode  at  x  =  0 
and  1. 


The  optimum  electron  beam  radius,  and  the  peak  gain  at  that  radius,  depend  strongly 
on  the  beam  emittance.  With  better  emittance,  curve  (a)  in  Figure  21,  the  optimum 
electron  beam  radius  is  small,  cse  -  0.3,  and  peak  gain  is  increased  to  G'G0  =  66%. 
Increasing  emittance  causes  the  optimum  oe  to  grow  larger  and  the  maximum  to 
become  more  broad  until  (curve  (c))  the  maximum  gain  is  only  weakly  dependent  on 
ce  due  to  the  large  angular  spread  at  any  radius.  This  maximum  gain  is  reduced  to 
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GIG o=  10%  at  cts  =  0.6.  The  expected  emittance  in  the  LPC  UVFEL  is  e  =  0.43, 
curve  (b).  We  find  that  the  optimum  electron  beam  radius  is  ce  =  0.4,  slightly  smaller 
than  that  proposed. 


Figure  21 .  Plots  of  normalized  gain  versus  electron  beam 
radius  for  three  different  beam  emittances:  (a)  e  =  0.28, 

(b)  e  =  0.43,  and  (c)  e  =  0.6. 

E.  CONCLUSIONS 

The  small  characteristic  mode  size  of  the  short  wavelength  UVFEL  makes  the 
electron  beam/optical  mode  overlap  a  crucial  consideration.  Making  the  optical  mode 
larger  by  either  decreasing  or  increasing  the  Rayleigh  length  from  z0  =  0.3  results  in 
increased  normalized  gain  for  a  given  electron  beam  radius.  This  is  an  important 
result  for  UVFELs.  It  shows  that  it  is  not  necessary  to  have  the  electron  beam 
enveloped  by  the  optical  mode  throughout  their  interaction.  With  Rayleigh  length  fixed, 
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focussing  the  electron  beam  to  a  smaller  radius  will  increase  gain  until  the  angular 
spread,  which  must  increase  as  the  electron  beam  radius  is  decreased,  degrades  the 
interaction  and  gain  is  reduced.  With  improvements  in  emittance,  the  optimum 
electron  beam  radius  is  smaller  and  the  peak  FEL  steady-state  gain  is  significantly 
increased.  We  find  the  optimum  electron  beam  radius  for  the  LPC  UVFEL  is  ce  =  0.4. 
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V.  FELS  WITH  SHORT  RAYLEIGH  LENGTH 


A.  MOTIVATION 

The  problem  of  high  power  density  on  optical  cavity  mirrors  is  shared  by  both 
UVFELs  and  high  average  power  FELs,  the  former  due  to  its  short  wavelength  and 
therefore  a  small  diffraction  rate,  the  latter  due  to  extremely  high  (=  100  MW) 
intracavity  average  power.  This  problem  could  be  solved  by  designing  a  long  optical 
cavity  with  z0  =  0.3,  allowing  space  for  diffraction  to  make  the  beam  spot  on  the 
mirrors  large.  Another  way  to  make  a  large  spot  is  to  use  a  cavity  with  a  short 
(z0  <  0.3)  Rayleigh  length.  When  the  FEL  has  limited  space  for  the  optical  cavity  (on  a 
ship,  for  example),  the  second  method  of  mirror  protection  is  particularly  attractive. 

As  an  example  of  the  effect  of  the  Rayleigh  length  on  mirror  power  density, 
consider  a  1  MW  FEL  at  X  =  1  pm  and  0  =  5.  Simple  calculations  based  on  (60)-(63) 
show  that  the  mirror  power  density  in  a  20  m  long  cavity  with  z0  =  0.2  is  about 
1500  kW /cm2,  which  is  quite  large.  However,  if  we  reduce  the  Rayleigh  length  to 
z0  =  0.02  and  shorten  the  cavity  length  to  1 0  m,  the  mirror  power  density  is  reduced  to 
about  600  kW/cm2.  That  is,  we  can  reduce  the  power  density  on  the  mirrors  by  more 
than  a  factor  of  two  and  cut  the  cavity  length  in  half  by  reducing  the  Rayleigh  length  by 
a  factor  of  ten. 

Reducing  the  Rayleigh  length  is  accomplished  by  making  the  mirrors  larger  and 
their  radii  of  curvature  smaller  (recall  z0  =  (2 SR  -  SZ)'A/2L),  allowing  the  light  to  cover 
a  large  mirror  while  focussed  to  a  small  waist  at  the  center  of  the  cavity.  The  stability 
criteria  for  a  symmetric  optical  cavity,  0  <  (1  -  S/R)(  1  -  S/R)  <1,  reduces  to  z|  >  0 
(Ref.  7).  Therefore,  as  long  as  there  is  a  Rayleigh  length  (i.e.,  z0  >  0),  the  optical 
cavity  will  be  stable. 

In  section  G  of  Chapter  II,  we  discussed  the  physics  of  the  FEL  interaction  in 
the  Gaussian  optical  mode.  We  discussed  how  the  optical  field  amplitude  laGl  is 
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greatest  at  i  =  xw,  which  causes  the  phase  space  separatrix  height  to  increase  as 
t  :  0  -»  xw  then  decrease  as  x  :  xw  -»  1.  This  effect,  which  is  proportional  to  i ir1/2(x), 
complicates  the  phase  space  paths  and  will  be  exacerbated  by  a  small  Rayleigh 
length.  The  short  Rayleigh  length  implies  a  small  mode  waist  where  the  optical  field 
amplitude  will  be  much  larger  than  at  the  ends  of  the  undulator.  We  also  saw  how  the 
phase  of  the  Gaussian  mode  caused  the  separatrix  to  shift  to  the  right  as  x  :  0  -»  1 . 
This  requires  the  initial  phase  velocity  for  maximum  gain  to  be  increased  from  v0  =  2.6, 
the  value  found  when  using  the  plane  wave  approximation.  For  z0  =  0.1,  the  phase 
shift  for  on-axis  electrons  is  A<|>  =  <|>(t=1  )  -  <|>(x=0)  =  0.9rc,  and  is  smaller  for  off-axis 
electrons.  With  a  small  Rayleigh  length,  the  wavefront  curvature  has  to  be  larger  in 
order  that  the  light  be  focussed  to  a  small  waist  and  rapidly  expand  to  cover  the 
mirror. 

Operating  the  FEL  with  a  short  Rayleigh  length  is  not  the  optimum  design.  The 
previous  paragraph  showed  how  the  interaction  of  the  electrons  and  the  light  changes 
with  Rayleigh  length.  This  concept  has  never  been  studied  before  now,  precisely 
because  the  design  is  not  optimum.  In  this  chapter  we  will  show  the  consequences  of 
using  a  short  Rayleigh  length.  Since  we  will  analyze  the  case  z0<1,  we  cannot 
expand  (60)  in  1/(z0)  as  we  did  in  Chapter  II,  but  must  use  simulation  instead.  We  will 
analyze  the  UVFEL  of  Chapter  IV  first,  both  the  single  pass  gain  G(x)  and  the  gain 
evolution  G(n)  over  several  passes  when  in  steady  state.  Next  we  present  a  MW- 
class  FEL  design  based  on  the  short  Rayleigh  length  concept  and  simulations  of 
same. 

B.  UVFEL  WITH  SHORT  RAYLEIGH  LENGTH 

We  first  test  the  short  Rayleigh  length  FEL  concept  by  simulating  the  LPC 
UVFEL,  which  has  parameters  listed  in  Table  2. 
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1.  Single  Pass  Gain 

We  study  the  single  pass  gain  G( t)  in  the  UVFEL  using  transverse  multimode 
simulations  like  those  discussed  in  the  previous  chapter.  In  this  case,  however,  since 
only  one  pass  is  simulated,  cavity  mirrors  do  not  enter.  The  optical  field  is  initialized 
as  a  Gaussian  mode  and  can  only  evolve  by  mode  distortion  caused  by  the  electron 
beam  within  the  single  pass.  Contrast  this  with  the  multiple  pass  simulations 
presented  in  the  previous  chapter  in  which  the  optical  field,  entrained  by  the  mirrors, 
evolves  self-consistently  over  many  passes. 

Plotted  in  Figure  22  is  the  single  pass  gain  evolution  G(t)  along  the  undulator 
for  four  different  electron  beam  radii:  (a)  oe  =  0.2,  (b)  ae  =  0.45,  (c)  oe  =  0.9,  and  (d) 
oe  =  1.6. 


T 


Figure  22.  Gain  along  the  undulator  G(t)  for  electron 
beams  with  radius  (a)  ae  =  0.2,  (b)  oe  =  0.45,  (c) 
oe  =  0.9,  and  (d)  ae  =  1.6.  Rayleigh  length  is  z0  =  0.3, 
mode  waist  radius  is  w0  =  0.5. 
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The  Rayleigh  length  is  z0  =  0.3  giving  a  mode  waist  radius  of  w0  =  0.5  (see  (63)).  The 
expected  LPC  UVFEL  electron  beam  emittance  e  =  0.43  is  again.  Electron  beams  (a) 
and  (b)  have  radii  smaller  than  the  mode  waist  radius  (i.e.,  oe  <  w0)  so  that  they  are 
enveloped  by  the  light  for  the  duration  of  their  interaction.  The  radii  of  electron  beams 
(c)  and  (d)  are  greater  than  the  mode  waist  radius.  Because  of  this,  their  final  gain 
G(t=1)  is  reduced  compared  with  beams  (a)  and  (b).  Also  note  that  because  of  the 
large  angular  spread  o0  =  4.6  for  electron  beam  (a),  the  final  gain  is  smaller  than  that 
for  beam  (b).  We  saw  in  section  D  of  the  last  chapter  that  electron  beam  (b)  with 
oe  =  0.45  and  oe  =  0.9  was  nearly  optimum  for  z0  =  0.5  (optimum  was  oe  =  0.4). 

The  most  encouraging  aspect  of  the  plots  in  Figure  22  is  that  although  the  final 
gain  differs  among  the  four  electron  beams,  all  show  nearly  exponential  gain  (see 
(58))  throughout  the  length  of  the  undulator.  This  point  is  better  illuminated  by  Figure 
23.  In  this  figure  we  normalized  the  curves  G(x)  from  Figure  22  to  their  respective  final 
gains  G(x=1)  in  order  to  compare  the  structure  of  the  curves  more  closely.  Figure  23 
shows  convincingly  that  all  four  of  the  electron  beams  simulated  resulted  in  similar 
gain  evolution  throughout  the  length  of  the  undulator,  even  those  ((c)  and  (d))  with 
oe  >  wQ.  We  see  that  the  gain  evolution  starts  in  all  four  beams  after  a  small  delay  of 
t  =  0.2,  which  is  just  the  time  it  takes  for  the  electron  bunch  to  form  (Ref.  3). 

In  Figure  24  we  have  plotted  G(x)  for  the  same  four  electron  beam  radii,  but 
with  the  Rayleigh  length  reduced  to  z0  =  0.1.  This  means  w0  =  0.3  and  only  electron 
beam  (a)  with  ae  =  0.2  fits  inside  the  optical  mode  throughout  their  interaction.  For 
this  reason,  beam  (a)  has  the  highest  final  gain.  Electron  beam  (b)  has  final  gain 
reduced  by  nearly  a  factor  of  two  over  that  in  Figure  22  because  it  no  longer  fits  inside 
the  optical  mode  at  the  waist.  Beams  (c)  and  (d)  also  have  final  gain  reduced,  but 
only  beam  (d)  has  final  gain  less  than  10%.  This  figure  shows  that  in  a  single  pass, 
where  the  electron  beam  cannot  couple  to  higher  order  modes,  it  is  more  important  to 
have  the  electron  beam  enveloped  by  the  light  than  to  have  a  small  angular  spread. 
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Figure  23.  Normalized  gain  along  the  undulator 
G(t)/G(t=1)  for  electron  beams  with  radii  (a)  oe  =  0.2,  (b) 
ce  =  0.45,  (c)  oe  =  0.9,  and  (d)  oe  =  1.6.  Rayleigh 
length  and  mode  waist  radius  same  as  in  Figure  22. 


For  Figure  25,  the  Rayleigh  length  was  further  shortened  to  z0  =  0.02.  In  this 
case  w0  =  0.14  and  none  of  the  four  electron  beams  fit  inside  the  optical  mode  at  the 
waist.  The  final  gain  for  electron  beams  (a)  and  (b)  is  less  than  half  their  value  at 
Zq  =  0.1  in  Figure  24.  However,  there  is  still  gain  at  this  very  short  Rayleigh  length, 
provided  the  electron  beam  is  small  enough.  We  see  that  oe  =  0.9  or  1.6,  curves  (c) 
and  (d),  are  clearly  too  large.  The  final  gain  G(t)  >  20%  for  the  two  smallest  electron 
beams. 
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Figure  24.  Gain  along  the  undulator  G(x)  for  electron 
beams  with  radius  (a)  oe  =  0.2,  (b)  oe  =  0.45,  (c) 
oe  =  0.9,  and  (d)  oe  =  1.6.  Rayleigh  length  is  z0  =  0.1, 
mode  waist  radius  is  w0  =  0.3. 


In  order  to  better  compare  the  evolution  of  gain  G(t)  in  Figure  25,  the  curves 
are  normalized  to  their  respective  final  values,  as  in  Figure  23.  The  result  is  plotted  in 
Figure  26.  The  effect  of  the  short  Rayleigh  length  is  clear:  single  pass  gain  increases 
greatest  where  the  fields  la  I  are  strongest,  at  the  beam  waist  x  =  0.5.  The  short 
Rayleigh  length  causes  the  fields  at  x  =  0.5  to  be  much  stronger  than  at  any  other 
position.  This  can  be  seen  in  electron  beams  (a),  with  ae  =  0.2,  and  (d),  with 
oe  =  1 .6.  Electron  beam  (a)  is  focussed  at  x  =  0.5;  it  was  designed  to  take  advantage 
of  the  field  strength  at  the  center  by  concentrating  electron  beam  current  there. 
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Electron  beam  (d)  has  a  radius  eight  times  that  of  beam  (a).  It  is  so  large  that  gain  is 
poor  throughout  the  interaction  at  this  Rayleigh  length.  What  little  final  gain  exists, 
however,  was  enhanced  at  t  =  0.5  where  optical  fields  are  strongest.  Indeed,  we  see 
in  Figure  26  that  there  was  net  absorption  for  beam  (d)  until  the  field  strength  started 
to  increase  at  i  ~  0.4.  Also  evident  in  the  figure  is  that  electron  beams  (b)  and  (c) 
have  radii  which  allow  gain  to  be  distributed  more  uniformly.  That  is,  they  are  larger 
than  beam  (a)  so  gain  does  not  increase  as  much  in  the  center  of  the  undulator,  but 
are  smaller  than  beam  (d)  so  that  gain  increases  more  at  the  ends  of  the  undulator. 


Figure  25.  Gain  along  the  undulator  G(x)  for  electron 
beams  with  radius  (a)  oe  =  0.2,  (b)  ae  =  0.45,  (c) 
a0  =  0.9,  and  (d)  ae  =  1 .6.  Rayleigh  length  is  z0  =  0.02. 
mode  waist  radius  is  w0  =  0.14. 
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Figure  26.  Normalized  gain  along  the  undulator 
G(x)/G( t  =  1)  for  electron  beams  with  radius  (a)  ae  =  0.2, 
(b)  ce  =  0.45,  (c)  ae  =  0.9,  and  (d)  oe  =  1.6.  Rayleigh 
length  and  mode  waist  radius  same  as  in  Figure  25. 


To  summarize,  we  have  found  that  using  a  very  short  Rayleigh  length  degrades 
gain.  However,  gain  is  still  significant  for  the  LPC  UVFEL  at  a  Rayleigh  length  as 
short  as  z0  =  0.02.  For  ae  =  0.2,  the  final  gain  is  G  =  37%,  while  for  ae  =  0.45,  the 
final  gain  is  G  =  23%.  The  gain  is  only  seriously  degraded  for  electron  beams  with 
oe  =  0.9  and  oe  =  1.6,  which  have  radii  two  and  four  times  greater  than  optimum, 
respectively.  We  also  found  that  the  most  important  region  for  the  interaction  is  the 
center  of  the  optical  mode  where  the  fields  are  strongest,  even  if  the  electron  beam  is 
much  larger  than  the  mode  waist.  This  makes  sense  because  it  is  the  field  amplitude 
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I  a  I  which  drives  bunching  in  the  pendulum  equation  (36). 

When  this  research  was  started,  however,  it  was  thought  that  the  interaction 
may  "turn  off"  at  the  center  of  the  undulator,  meaning  that  section  of  undulator  could 
be  replaced  by  a  drift  space,  making  an  optical  klystron,  which  would  improve  gain. 
The  optical  klystron  consists  of  two  undulator  sections  with  a  drift  space  at  the  center 
(Ref.  3).  The  first  section,  called  the  modulator,  starts  the  electron  bunching  process. 
The  drift  space  lets  the  electron  momenta  continue  bunching  in  the  absense  of  the 
undulator  field.  The  second  undulator  section,  called  the  radiator,  causes  the  bunched 
electrons  to  radiate  coherently.  Instead  of  a  drift  space,  the  same  effect  can  be 
obtained  in  a  shorter  distance  with  a  dispersive  section,  which  is  just  a  magnet  that 
bends  the  electron  beam  out  of  and  back  into  the  beam  path.  Slower  electrons  bend 
through  a  smaller  radius  than  the  more  energetic  electrons.  Since  the  faster  electrons 
travel  a  greater  distance  through  the  dispersive  section  than  the  slower  electrons,  they 
are  bunched  upon  exiting.  It  was  found,  however,  that  we  do  not  need  to  resort  to 
exotic  klystron  undulator  designs  to  get  the  desired  gain. 

2.  Steady  State  Gain 

Having  studied  how  the  gain  G(t)  evolves  within  a  single  pass,  we  will  now 
simulate  several  passes  and  follow  the  gain  evolution  G(n).  Recall  that  G(n)  follows 
the  final  gain  G  at  t  =  1  after  each  pass  n.  We  use  the  transverse  multimode, 
multiple  pass  simulations  used  in  the  previous  chapter  (see,  for  example,  Figure  20) 
where  the  optical  mode  evolves  self-consistently  in  the  optical  cavity  over  many 
passes. 

The  electron  beams  (a),  (b),  and  (c)  from  the  previous  figures  were  simulated. 
The  simulations  start  from  weak  fields  a0  =  0.01  and  evolve  until  the  limit  of  weak 
fields  la  I  =  it  is  reached.  At  this  point  the  final  gain  is  in  steady-state,  by  which  we 
mean  that  single  pass  gain  G(t=1)  has  nearly  the  same  value  over  several  passes 
and  G(n)  is  nearly  constant.  It  is  this  steady-state  value  of  the  gain  evolution  G(n) 
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that  is  plotted  versus  the  Rayleigh  length  z0  in  Figure  27. 


Figure  27.  Steady-state  gain  G  versus  Rayleigh  length 
for  electron  beams  with  radius  (a)  oe  =  0.2,  (b)  oe  =  0.45, 
and  (c)  oe  =  0.9.  We  find  G  >  60%  at  z0  =  0.01  for 
beams  (a)  and  (b). 


In  Figure  27,  we  see  that  the  nearly  optimum  electron  beam,  (b),  has  the 
highest  steady-state  gain  at  all  Rayleigh  lengths  explored.  This  is  in  contrast  with  the 
single  pass  studies  in  which  electron  beam  (a)  had  highest  final  gain  for  zQ  <  0.3. 
This  electron  beam  is  evidently  more  efficient  in  coupling  to  higher  order  modes  than 
any  of  the  other  beams  explored.  Perhaps  this  is  because  of  the  better  balance 
between  it  s  beam  radius  ce  =  0.45  and  angular  spread  oe  =  0.9  in  their  contribution  to 
the  electron  phase  velocity  spread  (i.e.,  it  is  better  "matched"  (Ref.  3)).  The  maximum 
gain  for  all  three  of  the  beams  occurs  at  Rayleigh  length  z0  =  0.1.  The  most  important 
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feature  to  note  in  the  figure  shows  the  goal  of  this  chapter:  steady-state  gain  is 
substantial  even  at  z0  =  0.01.  For  electron  beams  (a)  and  (b)  we  have  G(n)  >  60%, 
while  that  for  beam  (c),  twice  the  optimum  radius,  is  still  a  respectable  12%. 

Our  conclusion  is  that  the  FEL  with  a  short  Rayleigh  length  may  be  able  to 
operate  above  threshold  with  respectable  gain  while  spreading  the  optical  mode  at  the 
mirrors.  The  trade-off  is  slightly  lower  gain  for  a  much  shorter  optical  cavity  with 
reduced  power  loading  at  the  mirrors. 

C.  MW-CLASS  FEL  WITH  SHORT  RAYLEIGH  LENGTH 

Encouraged  by  success  in  modelling  the  UVFEL  with  a  short  Rayleigh  length, 
the  concept  was  used  to  design  a  MW-class  FEL. 


1.  MW-Class  FEL  Design 

Our  goal  in  designing  a  MW-class  FEL  is  to  derive  a  set  of  parameters  which 
will  provide  sufficient  gain  within  an  optical  cavity  of  length  10  m.  This  length  was 
chosen  because  it  fits  within  the  beam  of  a  modern  naval  cruiser  or  destroyer. 

FEL  gain,  discussed  in  section  E  of  Chapter  II,  is  dependent  upon  the 
dimensionless  current  density  j  (see  (57)  and  (58)).  Maximizing  gain  therefore 
requires  that  we  maximize  j.  Recall  that  the  definition  of  the  current  density  in  the 
linear  undulator  is 


.  6N(enK(J0^)-J^))Lfp 

1  yW 

where  £  =  K^/(2(1  +  K2)).  Since  N  =  L/\q  we  can  re-write  (69)  as 
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where  we  have  separated  the  parts  of  j  according  to  their  source.  Following  the 
fundamental  constants  is  the  undulator  contribution  K2(JJ)2L?I\0  where  we  have 
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shortened  the  Bessel  function  term  in  (70)  to  JJ.  Since  K  ~  1  in  an  FEL,  which  also 
fixes  the  Bessel  terms,  we  are  left  with  L  and  X0  as  designable  undulator  parameters 

(Ref.  3).  The  electron  beam  contribution  is  p/y3.  The  charge  density  p  depends  on  the 
peak  current  t  through 


Jiecr2(/e/c)  '  ' 

where  Q  =  t  x  (/e/c)  is  the  charge  per  micropulse  (not  the  cavity  quality  factor),  re  is 
the  electon  beam  radius,  and  le  is  the  length  of  the  micropulse.  The  electron  beam 
energy  y me2  determines  y.  Therefore  we  can  write 
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Recall  that  the  emission  wavelength  is  given  by 


A.  —  A,0 


(1  +  *2) 
2Y2 


(73) 


With  K  -  1  this  becomes  X  =  Vt2-  This  means  that  if  we  fix  the  wavelength  X,  then 
X0  and  y  are  not  independent.  Therefore,  the  remaining  free  parameters  are 
^  i  A-q,  and  X,. 


Our  philosophy  in  choosing  a  design  was  to  use  a  short  undulator  and  use 
peak  current  to  maximize  j.  The  reasons  for  this  are  twofold.  First,  we  are  using  an 
optical  cavity  10  m  long  with  a  short  Rayleigh  length.  The  undulator  has  to  be  short  in 
order  to  fit  around  the  optical  mode  at  the  center  yet  not  "scrape"  any  of  the  light, 
which  would  damage  the  magnets.  Second,  a  short  undulator  means  the  FEL  will  be 
less  susceptible  to  electron  beam  quality.  This  is  because  the  gain  spectrum 
bandwidth,  which  is  proportional  to  1/(2 N),  increases  for  decreasing  undulator  length 
L  =  NXq.  This  was  discussed  in  section  E  of  Chapter  II  and  is  important  to  the  design, 

because  as  we  demand  more  current  from  the  accelerator  we  expect  the  beam  quality 
to  suffer. 
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Accordingly,  we  chose  an  undulator  length  L  =  1m.  The  undulator  wavelength 
was  taken  to  be  X0  =  4  cm,  which  means  N  =  25  periods.  Since  the  wavelength 
X  =  1  pm  is  fixed  by  the  ship  self-defense  application,  the  electron  beam  energy  is 
ymc2  =  ^X0/Xmc2  =100  MeV.  For  the  remaining  parameters  we  chose  f  =  600  A, 
re  =  0.3  mm,  and  le/c  =  3  ps.  This  results  in  a  total  charge  per  micropulse 
Q  =  t  x  ( le/c )  =  1.8  nC.  If  every  RF  bucket  of  an  f  =  500  MHz  accelerating  field  is 
filled  by  the  electron  injector,  the  duty  factor  is  (/e/c)  x  f  =  0.0015.  The  average 
current  is  then  7=  0.0015  x /  =  0.9  A.  This  results  in  an  average  electron  beam 
power  of  Fe  =  TV  =  0.9  A  x  IOOxIO6  V  =  90  MW.  An  efficiency  of  1.1%  will  be 
needed  to  extract  1  MW  from  this  beam.  These  design  parameters  resulted  in  y  =  12. 
All  parameters  are  summarized  in  Table  3. 

For  electron  beam  quality  we  anticipate  a  trade-off  between  high  peak  current 
and  electron  beam  quality.  We  assumed  an  energy  spread  Ay/y  =  0.5%  and  a 
normalized  emittance  ew  =  20n  mm-mrad.  These  choices  resulted  in  cq  =  1 .5  and 
e  =  0.27.  The  dimensionless  electron  beam  radius  is  ae  =  0.5  so  that  the  angular 
spread  must  be  oe  =  0.3.  Noting  that  N  is  included  in  the  numerator  of  (18)  and  (56), 
we  see  that  the  small  values  for  oe ,  ae,  and  e  reflect  our  choice  of  a  short  undulator. 
The  beam  quality  parameters  are  also  summarized  in  Table  3. 

2.  Simulations  of  the  MW-Ciass  FEL 

Now  that  we  have  derived  our  parameters  for  the  MW-class  FEL,  we  use 
simulations  to  make  sure  it  works.  Transverse  multimode,  multiple  pass  simulations 
using  the  parameters  in  Table  3  were  used  to  obtain  the  steady-state  gain  at  different 
Rayleigh  lengths.  An  example  of  these  simulations  at  z0  =  0.1  is  shown  in  Figure  28. 
Steady  state  gain  is  G(n)  =  26%,  plotted  at  lower  right,  and  power  P(n)  is  increasing 
approximately  exponentially.  The  optical  field  evolves  over  20  passes  yet  is  still  weak 
(la  l  c  n).  Bunching  is  not  evident  in  the  phase  space  plot  at  lower  center. 
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Beam  Energy 

100  MeV 

Peak  Current 

600  A 

Average  Current  (@500  MHz) 

0.9  A 

Micropulse  Length 

3.0  ps 

Charge  per  Micropulse 

1.8  nC 

Electron  Beam  Radius 

0.3  mm 

£w 

20ji  mm-mrad 

Ay/y 

0.5% 

Undulator  Type 

Linear 

K 

1.0 

N 

25 

*0 

4.0  cm 

X 

1 .03  pm 

i 

12 

0.5 

°e 

0.3 

e 

0.27 

1.5 

Table  3.  Design  parameters  for  a  MW-class  FEL  suitable 
for  ship  self  defense. 
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FEL  WAVEFRONTS 
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Figure  28.  Transverse  multimode  simulation  of  MW-class 
FEL.  Gain  evolution  G(n)  is  in  steady  state. 


The  electron  positions  overlaid  on  the  plot  of  the  optical  field  amplitude  between  the 
mirrors  la (x,t)I  show  that  the  electron  beam  does  not  expand  nor  focus  as  x  :  0  -»  1. 
This  is  because  the  angular  spread  o6  =  0.3  is  relatively  small. 

Many  simulations  such  as  that  in  Figure  28  were  used  to  find  the  effect  of  the 
Rayleigh  length  on  steady-state  gain.  The  result  is  plotted  in  Figure  29.  As  in  our 
previous  analysis  of  the  UVFEL,  maximum  gain  G  ~  23%  for  this  FEL  occurs  at 
z0  =  0  1  (  see  Figure  28).  Even  at  the  short  Rayleigh  length  z0  =  0.01,  the  steady- 
state  gain  G  =  18%  is  sufficiently  large  that  simulations  indicate  the  design  will  work. 
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Figure  29.  Steady-state  gain  G  versus  Rayleigh  length 
z0  for  MW-class  FEL.  Peak  gain  G  =  23%  occurs  at 
z0  =  0.1.  At  z0  =  0.01,  gain  is  G  =  18%. 


D.  CONCLUSIONS 

Using  the  LPC  UVFEL  parameters  we  showed  that  the  concept  of  the  short 
Rayleigh  length  FEL  was  viable.  This  had  never  been  researched  before.  Based  on 
those  results,  parameters  for  a  MW-class  FEL  for  ship  self-defense  were  derived.  This 
FEL  design  was  simulated  and  steady-state  gain  was  found  to  be  G  =  18%  at 
Zq  =  0.01.  It  was  presented  at  the  workshop  entitled  "Navy  MW-Class  SSD  FEL 
Concepts  at  TJNAF  in  Newport  News,  Virginia,  which  was  sponsored  by  the  Navy 
High  Energy  Laser  office  at  SPAWAR.  It  is  one  of  only  two  designs  that  survived  the 
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workshop. 

To  improve  on  our  design,  the  electron  beam  radius  should  be  optimized.  At 
oe  =  0.5,  we  have  a  small  oe  =  0.3  which  suggests  it  would  be  advantageous  to  focus 
the  electron  beam  further.  Research  should  also  determine  what  affect  the  Rayleigh 
length  has  on  FEL  efficiency.  At  short  Rayleigh  lengths,  many  electrons  are  outside 
the  optical  mode  at  the  waist  and  do  not  participate  in  the  interaction.  Efficiency  might 
be  degraded.  This  is  the  trade-off  made  in  choosing  a  short  Rayleigh  length  over 
more  traditional  FEL  designs.  As  we  noted  above,  we  must  extract  1.1%  of  the 
electron  beam  energy  in  order  to  obtain  1  MW  average  power.  Methods  of  improving 
efficiency  should  also  be  studied. 

The  design  philosophy  of  using  a  shorter  optical  cavity  and  undulator  requires 
high  peak  and  average  current.  It  shifts  the  technology  burden  from  the  mirrors  to  the 
accelerator,  where  improvements  in  high  average  current  generation  and  acceleration 
must  be  made.  Another  area  of  accelerator  physics  that  needs  much  further  research 
is  in  electron  beam  transport.  For  our  laser  to  work  on  a  ship,  we  have  determined 
that  we  must  use  an  energy  recovery  architecture  similar  to  that  proposed  for  the  LPC 
UVFEL  in  order  to  increase  overall  system  efficiency,  among  other  reasons  (Ref.  9). 
This  reduces  the  power  consumption  of  the  FEL  by  recirculating  the  electron  beam 
through  the  accelerator,  where  it  is  decelerated,  giving  up  its  energy  to  the  RF 
accelerator  field.  Recirculation  requires  the  bending  of  a  high  current,  high  energy 
electron  beam  through  two  180  degree  bends,  but,  at  the  time  of  this  writing,  there  is 
some  question  as  to  whether  or  not  that  is  even  possible. 
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VI.  COMPTON  BACKSCATTERING  OF  LASER  LIGHT 


A.  PRELIMINARIES 

We  turn  our  attention  now  to  the  scattering  of  laser  light  from  a  high  energy 
electron  beam.  The  electron  beam  scatters  high  energy  photons  out  of  an  incoming 
laser  field.  The  process  is  similar  to  the  FEL  in  that  the  incoming  laser  pulse  can  be 
thought  of  as  an  electromagnetic  undulator,  but  with  extremely  small  wavelength.  The 
process  has  been  described  classically  and  with  relativistic  quantum  mechanics, 
assuming  a  plane  wave  approximation  for  the  laser  field  (Refs.  11,  12,  13).  Recently, 
a  correction  to  the  plane  wave  approximation  was  attempted,  using  quantum 
electrodynamics  (QED)  (Ref.  14).  Derlet  et.  al.  treated  the  interaction  of  free  electrons 
with  rays  of  laser  photons  focused  by  a  lens  to  the  interaction  region  (Ref.  14).  The 
calculation  presented  in  this  chapter  is  the  first  attempt  to  correct  the  plane  wave 
approximation  to  the  laser  beam  by  treating  it  as  a  diffracting  beam.  That  is,  we 
assume  the  incoming  laser  beam  is  a  fundamental  Gaussian  mode. 

In  the  next  section,  we  calculate  the  vector  potential  for  the  Gaussian  optical 
mode.  We  approximate  the  Gaussian  mode  for  long  Rayleigh  length  in  a  manner 
similar  to  the  analysis  in  section  G  of  Chapter  II.  In  section  C,  we  calculate  electron 
wavefunctions  in  the  field  represented  by  our  vector  potential.  We  use  those 
wavefunctions  in  section  D  to  find  the  distribution  and  rate  of  photon  emission.  This 
will  be  the  first  time  this  calculation  has  been  done  without  the  plane  wave 
approximation  to  the  laser  field. 

The  notation  used  in  this  chapter  is  that  of  Sakurai  (Ref.  15).  We  represent  a 
four-vector  as,  for  example,  bM  with  p  =  1,2, 3.4  where 

bM  =  (b„b2,b3,b4)  =  (B,ib0)  (74) 

where  £>,.  b2.  and  b3  are  real,  corresponding  to  the  x-,  y-,  and  2-components  of  the 
three-vector  B,  and  b4  -  ib0  is  imaginary.  We  use  a  summation  convention  so  that 
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(75) 


the  scalar  product  of  two  four  vectors  b^  and  c ^  is 

3 

^  ~  b —  bftCQ  . 

/= i 

We  will  frequently  shorten  the  notation  for  a  scalar  product  from  b^c^  to  b  e  for 
brevity.  In  our  summation  convention,  Greek  indices  such  as  p,  v,  etc.  run  from  1  to  4 
and  Roman  indices  like  /,  j,  k,  etc.  run  from  1  to  3.  We  define  the  coordinate  vector 
as 

*n  =  C >t,ict )  ,  (76) 

where  It  =  x^X  +  x2$  +  x32.  The  momentum  four-vector  is 

p„  =  ( p,iE/c )  ,  (77) 

where  p  and  E  are  the  particle  momentum  and  total  energy,  respectively.  The 
wavevector  of  a  photon  is 

k n  =  (fr  ,/co/c)  ,  (78) 

where  Ik  I  =  2tz/X  =  co/c  for  a  free  photon.  Variables  are  not  made  dimensionless  as  in 
the  previous  chapters. 

We  use  Heaviside-Lorentz  units  throughout  this  chapter.  The  fields  and 
potentials  in  this  system  are  related  to  those  in  Gaussian  units  by  the  factor 
(Ref.  15).  We  will  take  the  electron  charge  e  to  be  negative  e  =  -le  I. 

A  typical  experiment,  and  the  one  we  will  use  in  order  to  obtain  order-of- 
magnitude  estimates  to  certain  terms,  is  the  proposed  "laser  synchrotron  source" 
(LSS)  at  the  Naval  Research  Laboratory  (NRL)  (Ref.  11).  They  have  proposed  using 
a  41  MeV  electron  beam  with  micropulse  length  le/c  =  1  ps.  The  incident  laser  is  a 
tabletop-terawatt  (T3)  device  with  wavelength  X  =  1  pm.  The  laser  peak  power  is  10 
TW  in  a  2  ps  pulse,  which  means  the  energy  per  pulse  is  20  J.  The  Rayleigh  length  is 
*o  =  7.9  mm  and  the  mode  waist  radius  is  w0  =  50  pm.  The  energy  of  the  emitted 
photons  is  expected  to  be  30  keV  (x-rays).  We  can  then  calculate  some  typical  values 
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for  useful  combinations  of  parameters:  eA0/mc 2  =  0.4,  KHmc 2  =  0.06,  and 
tko/mc2  =  2x1 0-6.  We  define  A0  as  the  amplitude  of  the  vector  potential,  m  is  the 
electron  rest-mass,  c  is  the  speed  of  light  in  vacuum,  /?  is  Planck’s  constant  divided  by 
2n,  £2  is  the  frequency  of  the  emitted  photon,  and  ro  is  the  frequency  of  the  incident 
laser.  Note  that  the  combination  eA0/mc2  is  analogous  to  the  undulator  parameter  K 
in  the  FEL. 


B.  VECTOR  POTENTIAL  FOR  LASER  LIGHT 


The  fundamental  Gaussian  mode  for  some  freely  propagating  scalar  field  U 
can  be  written: 


U  = 


Ur 


w(z) 


e 


x2+y2 

yw. 


(79) 


where  U0  is  the  amplitude  in  the  center  of  the  beam  waist  (Ref.  7).  The  mode  radius 
is 


w2(z)  =  w§ 


1  + 


(80) 


where  w0  is  the  radius  of  the  mode  waist.  The  mode  waist  radius  and  Rayleigh  length 
are  related  by 


The  optical  phase  $  is 


*o  = 


nw§ 


+  k 


x2  +  y2 
2  R(z) 


(81) 


(82) 


where  R(z)  =  z  +  z2/z  is  the  wavefront  radius  of  curvature.  The  origin  z  =  0  is 
located  at  the  mode  waist. 
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Since  our  goal  is  to  obtain  a  correction  to  the  plane  wave  approximation  for  our 
scattering  problem,  we  will  find  the  form  of  (79)  when  it  nearly  represents  a  plane 
wave,  or  z0  »  z.  For  this  case,  it  makes  sense  to  expand  (79)  in  the  small  parameter 
(1/z0)  in  a  manner  similar  to  the  WKBJ  approximation  (Ref.  16).  That  is,  we  will 
expand  both  the  amplitude  and  phase  of  (79)  in  (1/z0)  arid  approximate  them  by 
keeping  terms  only  to  first  order  in  (1/z0). 

Upon  expanding  the  amplitude  we  find  that 


U  o 


-e 


w(z) 

Likewise  we  find  for  the  phase 


x2+y2 

=  iJo’ 

.  w2(z). 

w0 

1 _ 1  /r(x2  +  y2) 

2 


+  O 


<|>  =  +  O 

zo 


Therefore  the  WKBJ  approximation  for  the  Gaussian  mode  is 


U  = 


Wo 


■\  _  1  k(x2  +  y2) 
zo  2 


-i[{k  -  Mz0)z-wt] 


(83) 


(84) 


(85) 


To  lowest  order  in  (1/Zq)  we  see  that  the  Gaussian  mode  is  nearly  a  plane  wave.  The 
only  differences  between  (85)  and  a  plane  wave  are  a  phase  shift  from  k  -*  k  -  1/z0 
and  an  amplitude  correction  which  decreases  the  amplitude  for  electrons  off-axis.  To 
an  electron  in  the  Gaussian  mode,  the  phase  shift  appears  as  a  shift  in  the  laser 
wavelength  X  -*  X/(1  -  1/kz0)  =  X(1  +  1//cz0).  Therefore  we  can  anticipate  that  the 
presence  of  the  Gaussian  mode  should  result  in  a  shift  in  wavelength  of  the  emitted 
radiation  as  well. 


For  our  calculation,  we  will  take  (85)  to  be  our  vector  potential.  Since  this 
potential  is  real  we  will  use  the  real  part  of  (85)  and  write 


A06 


1  ~  7^~^~2+  ^  cos[(fc  -  1/Zo)Z-Of) 


(86) 
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where  A0,  with  w0  absorbed,  is  the  on-axis  magnitude  of  A  and  a  is  the  polarization 
unit  vector.  If  the  laser  light  is  linearly  polarized  in  the  x-direction,  as  we  will  assume 
in  this  chapter,  then  3  =(1,0,0).  We  will  only  consider  on-axis  (i.e.,  x2  +  y2  =  0) 
electrons,  so  that  we  can  write  (86)  as  a  four-vector  as 

\(x)  =  (*(x),Mo)  =  ^Oa|icos(/Cvxv)  (87) 

where  =  (3,0)  and 

Kv  =  /fv  -  8V  =  0,0,k-—,i—  (88) 

z0  c 

is  the  modified  propagation  vector.  The  vector  8V  =  (0,0,1/z0,0)  represents  the  phase 
shift  discussed  above.  We  note  that  A^x)  satisfies  the  Lorentz  condition 
^A^/^x^  =  K^A^  =  0  and,  since  in  the  gauge  we  have  chosen  A0  =  0,  the  transversality 
condition  V-A  =  ft  /t  =  0  is  also  satisfied  (Ref.  15). 


C.  DIRAC  WAVEFUNCTIONS 


The  wavefunctions  for  an  electron  in  a  plane  wave  electromagnetic  field  were 
found  in  1935  by  D.  M.  Volkov  (Ref.  17).  We  will  follow  nearly  the  same  technique, 
but  using  our  vector  potential  (87)  for  the  Gaussian  mode. 

The  Dirac  equation  for  an  electron  in  an  electromagnetic  field  characterized  by 
A^ix)  is  given  by 


iiY4  +M 
»cYv  v  fi 


V  =  0 


(89) 


where  the  4x4  matrices  yv  are  called  the  "gamma  matrices"  (Ref.  15).  They  are 
defined  by  the  requirement  that  they  satisfy 

lYP.  Yv)  =  YnYv  +  YvYm  =  8Mv  (90) 

where  {  ,  )  is  called  the  anti-commutator,  and  8MV  =  1  if  p  =  v,  0  otherwise,  is  the 
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Kronecker  delta  (Ref.  15).  In  our  notation,  the  gamma  matrices  are 


0  -iok 


ia«  v  _  [  /  O' 

0  ’  w  -  [o  -I]  ’ 


where  /  is  the  2x2  identity  matrix  and  the  a*  are  the  Pauli  spin  matrices  given  by 

°'  =  [l  o]  '°2=[/  o\  -Os-to  -1.  • 

For  example,  we  write  y3  using  a3  as 


0  0  -/  O' 

V3=  ?  8  °0j,  .  (93) 

0  -/  0  of. 

We  also  note  here  that  the  gamma  matrices  are  Hermitian  yf=y4  (Ref.  15).  The 
wavefunction  in  (89)  is  a  four  component  object  called  a  Dirac  spinor  (Ref.  15).  Thus 
(89)  is  actually  four  coupled,  first-order,  differential  equations  for  \|/. 

We  can  obtain  the  quadratic  form  of  (89)  by  multiplying  on  the  left  by  the 
operator 

<94> 

with  the  result 


a 

a 

2/e.  a 

ie 

—  V 

a/Av' 

e2  .  . 

me 2 

axM 

/)c  MaxM ' 

hcYtJv 

.ajrpj 

tfc2  » 

n2  J 

xj/  =  0  . 


Substituting  (87)  for  Av  we  find  that 


- -  =  K  A  ’ 

axM 


where 


Av'  =  ->40av5in(Koxo) 


We  can  then  rewrite  (95)  as 
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a  a 


axH  ax^ 


- 'JL(y.K)(vA') _ ^->A2  -  mc2 

tic A*  ax4  t>c(y*)(yA}  rfc2  tf 


\|T  =  0 


(98) 


where  we  have  used  y-K  =  y^K^,  etc. 

Following  Volkov,  we  seek  a  solution  to  (98)  of  the  form 


iBr 

y  =  e  77  <!>(/c-x) 


(99) 


where  <|>(k-x)  is  a  Dirac  spinor.  We  subject  our  solution  to  the  initial  condition 
\|r  ->  u^(p)exp(ip-x/tl)  as  x  -*  -»  for  an  electron  coming  into  the  interaction  region. 
We  recognize  u{s)(p)  as  the  free  electron  spinor,  s  =  1 ,2  is  the  electron  spin  index 
and  p  is  the  electron  momentum  (Ref.  15).  The  initial  condition  means  physically  that 
prior  to  entering  the  interaction  region  the  electron  is  free.  Substituting  (99)  into  (98), 
we  find  that  §(k-x)  must  satisfy 


2/ 


PM'®  + 


^A-p-f^y-A')-^ 


m  =  o 


(100) 


where  i,  =  kx,  <t>'  =  a<})/a^,  and  we  have  used  k2  =  0.  This  first-order  differential 
equation  has  the  solution 


4>(S)  =  exp 


e 

2c 


j  CrJQ(r*0  d%. 


ie 

tic 


P-k 


e  A‘ 


2c  p-k 


u 


(101) 


where  u  is  a  Dirac  spinor.  Invoking  the  initial  condition,  we  find  that  <|>  -»  u  as  ^  ->  -« 
so  that  \j/ -*  uexp(ipxiti).  Therefore,  we  will  take  u  to  be  the  free  electron  spinor 
u  =  vmc^/EVu^ip)  where  V  is  the  normalization  volume  and  E  is  the  total  energy 
(Ref.  15).  The  wavefunction  for  the  incoming  electron  is 


V*.  =  exp 


2c  i. 


Ajj e  A _ 

p-k  2c  p-k 


(102) 


mc ‘ 


EV 


*  jEJL 

u{s\p)e  * 
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We  can  show  that  this  wavefunction  satisfies  the  Dirac  equation  (89),  as  it 
must,  in  addition  to  being  a  solution  to  (95)  (Ref.  12).  Let  D  represent  the  Dirac 
equation  operator  such  that  the  Dirac  equation  is  Dy  =  0.  Then  DD\ \t  =  0  represents 
(95).  (Recall  that  D  was  defined  in  (94).)  Let  %=Di|/;n  where  \]//n  is  given  by  (102). 
Now  since  \\ijn  reduces  to  the  free  electron  spinor,  which  satisfies  the  Dirac  equation, 
as  x  -co,  we  know  that  z->0asx^-».  Since  DDyjn  =  D%  =  0  and  %  ->  0  as 
*o  ->  we  must  have  %  =  0  identically.  Therefore,  Dyin  =  0  and  \\tin  satisfies  the 
Dirac  equation  (89). 

To  find  the  wavefunction  for  the  outgoing  electron,  we  must  solve  (95)  subject 
to  y  ->  u^\p')exp(ip'-x/Ti)  as  x0  -»  The  result  is 


Voir  =  exp 


oo 


ie  r  A  p' 
Tic  |  p'  ■ k 


e  <42 
2c  p'-k 


(103) 


x 


1/2  .pjx 

u^\p')e  * 


where  the  outgoing  electron  is  specified  by  spin  index  s',  momentum  p',  and  energy 
E . 

We  now  turn  to  the  task  of  evaluating  the  integrals  buried  inside  the 
wavefunctions  (102)  and  (103).  The  vector  potential  in  the  integrals  of  (102)  and  (103) 
is  a  function  of  K-x  =  k-x  -  8-x  while  the  variable  of  integration  is  k-x.  Let 
X  =  K  x  =  £-6x  and  perform  a  change  of  variables  to  %■  We  find 

r  1  r  (  ^  zr 

<^  =  1  +  ^T(6  x)  dx=  1  +  TT  +  0  T~  d *  (104) 

,  kz0  z0 

and  we  approximate  the  differential  d\  with  the  first  two  terms.  Making  the  change  of 
variables  in  (102)  we  write 
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V/n  =  exp 


(105) 


1  + 


1 

kz0 


exT{rK)(rA'(X'))dy. 

2 c  _i  p-k  X 


X+8'Xf 

je  f  AjzlP 

Kc  jL  [  P-k 


.?  A2(X') 

2c  p-k 

j 


1/2  fP-X 

u<-s\p)e  n 


The  correction  term  (1  +  Vkz0)  which  we  discussed  in  the  last  section  appears 
explicitly  in  the  electron  wavefunction.  Substituting  the  vector  potential  (87)  into  (105) 
and  performing  the  integrations  we  arrive  at 


V/n 


=  extf 


1  + 


e 

2c 


(rK)(ra) 

p-k 


cos  (k-x) 


(106) 


te 

f)c 


a-p 

p-k 


s\n(k-x) 


eA§ 

1 

-  -\N 

k-x  sin(2/cx) 

2c 

p-k 

<  J 

2  4  J 

%  .px 

u^s\p)e  n 

where  we  have  dropped  divergent  sine  terms  in  the  phase  of  the  wavefunction.  The 
divergence  of  these  terms  is  not  physical,  arising  because  our  vector  potential 
describes  a  laser  beam  undamped  at  while  we  require  that  the  fields  associated 
with  Ay  vanish  at  Likewise,  we  perform  the  change  of  variables  (104)  and  the 
integration  in  (103)  with  the  result 
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Vout  =  exn 


fi+  1 1 

e 

(y-K)(ya) 

[  kz°\ 

[2c 

p'-/c 

cos  (k-x) 


(107) 


le 

tic 


a-p' 


p’-k 


s\n(k-x)  - 


eA$ 


2c 


p'-k 


k-x  sin(2ft-x) 

2  4 


x 


mc‘ 


EV 


>/2 


u{s,)(p')e 


,p'-x 


for  the  outgoing  electron. 

We  now  have  the  wavefunctions  for  the  electron  entering  the  interaction  region 
(106)  and  for  the  electron  leaving  the  interaction  region  (107).  The  difference  between 
these  wavefunctions  and  the  Volkov  wavefunctions  is  in  the  factor  (1  +  1/(/rz0))  in  the 
exponential  and  the  presence  of  the  modified  propagation  vector  K instead  of  k M 
(Refs.  12,  13,  17).  As  we  expected,  by  choosing  a  Gaussian  mode  which  differed  only 
slightly  from  the  plane  wave  (i.e.,  z0  »  z),  our  wavefunctions  differ  only  slightly  from 
the  Volkov  wavefunctions.  We  now  use  these  wavefunctions  to  find  the  effect  of  the 
Gaussian  mode  on  photon  emission. 


D.  PHOTON  EMISSION 

For  the  calculation  of  the  differential  transition  rate  in  this  section  we  will  use 
"natural  units"  ti=  c  =  1.  This  aids  the  calculation  in  that  we  will  not  have  to  keep 
track  of  all  of  these  fundamental  constants.  We  will  write  our  results  at  the  end  of  the 
section  with  these  constants  included  so  that  real  numbers  may  be  calculated. 

The  amplitude  for  an  electron  transition  from  y(„  to  with  the  emission  of  a 
photon  is 

M  =  -\j»A»nt)dAx  (108) 
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where  the  integral  is  over  all  space  and  time  (Ref.  15).  The  interaction  vector  potential 
A^nt)  for  the  emission  of  a  single  photon  is 

AHnt)  =  V^7ef)e“/(/x)  (109) 

where  / ^  is  the  emitted  photon  propagation  vector  and  e^,  r  =  1,2,  is  its  polarization 
vector  (Ref.  15).  These  are  given  by 

V  =  (t  iO)  (110) 

where  Q  is  the  frequency  of  the  emitted  photon  and 

4r]  =  (4r)’°)  on) 

is  chosen  such  that  (£(1),&(2),7^)  form  a  right-handed  coordinate  system  in  three 
dimensional  space  (Ref.  15).  The  interaction  of  the  electrons  with  this  unquantized 
vector  potential  (109)  for  photon  emission  is  rigorously  equivalent  to  their  interaction 
with  the  quantized  electromagnetic  field  (Ref.  15).  The  transition  current  is 

=  out~f^  in  (112) 

where  is  given  by  (91)  and  y  =  v|/*y4  is  the  adjoint  spinor  (Ref.  15).  Substituting 
(109)  into  (108)  we  can  write  the  transition  amplitude  as 

M = ‘  (113) 

To  find  the  transition  current  yM(x)  we  substitute  the  wavefunctions  (106)  and 
the  adjoint  of  (107)  (recall  ft  =  c  =  1)  into  (112)  with  the  result 
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J '*(*>= Hw^ikel(p-'','x 


(114) 


x  l  v{*\p')exp  ^°i+-L  XMXti  cosk.x  y 
2  kz0  p'  k  ^ 


x  exp  1  +  rKr*  cos k-x  u{s)(p) 
2  kz0  p-k 


x  exp  msin/c-x  -  /'p(2/c-x  +  sin2/c-x) 


where  we  have  defined 


a  =  eAn  1  + 


1  a-p  _  a-p' 

/cz0  p  -/r  p'  k 


(115) 


e2/\  2  r  if  . 

8  kZn  D- 


kz0  p  ■ k  p'  k 


(116) 


We  will  now  simplify  and  make  some  approximations  to  (114).  Using  the 
property  of  gamma  matrices  that 


(yb)(yb)  =  b2  , 


(117) 


which  is  a  special  case  of 


(yb)(yc)  =  b  e  -  (yc)(yb)  , 


(118) 


we  can  expand  the  real  exponentials  in  (114).  The  expansion  yields 
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exp 


eA0 

r  i 

1  +  1 

(Y-a)  (yK)  cosk.zl 

2 

kz0 

p'-k 

(119) 


1  + 


eA( 


1  + 


1 

(y-a)(y-K) 

kz0 

p'-k 

where  we  have  ignored  terms  of  higher-order  in  (1/z0)  and  those  proportional  to 
( eA0/(2p')2(Vkz0 )  =  10-10  and  higher.  We  used  the  NRL  parameters  discussed  in 
section  A  to  calculate  these  estimates.  We  can  simplify  (114)  by  writing  the  complex 
exponential  involving  a  and  (3  in  terms  of  Bessel  functions  using  the  generating 
function 


e*sin e=  J  eim6Jm  (x)  (120) 

m-—o° 

where  m  is  an  integer  (Ref.  18).  We  find 

exp[/asin/c-x  -  /’Psin2fr-x]  =  £  £  e'>-2n>**Jm(a)Jn(p)  (121) 

m  =—oon  =— oo 


=  Z  el2qkx  Z  JZ(q+n)(o)J„M 

Cf  =  — oo  n  =— oo 

=  £  e/2‘**G0(g,a,p) 

Q=-oo 


where  G0(g,a,P)  =  £  ^2(q+n)(a)^n(P)-  This  term  (121)  leads  to  emission  into 

n=~<™ 

harmonics  where  q  is  the  harmonic  number.  We  then  incorporate  these  simplifications 
and  re-write  (114)  as 


x 


♦ e4° 

1  +  -1- 

J^cosk-j 

2 

kzD 

ff-k 

+ 

*  > 

^-ticoskx! 

2 

kz0 

V  > 

p-k 

Tm 

u(s>(^)  . 


(122) 
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Substituting  (122)  into  (113)  to  get  the  transition  amplitude  we  find  that  we  can 


write 


M  —  M0  +  M  i  +  M2 


(123) 


where 


x  U{*\p')rz{r)u{s\p)  G0(q,a,$) 


(124) 


and 


M- 


= -WisWlW  fv£jd'x  e^ptp’~2{q 


-W)x 


(125) 


eA( 


1  + 


kZr 


n(^\p') 


ya  y-K  ye(r)  ye^y/Cya 


P'Ar 


p-k 


u{s\p) 


x  G0(g,a,P)cos/r-x 


and  M2  «  (e/4o/m)2((o/m)  =  10~7  can  be  neglected  in  comparison  with  /W0.  For  the 
transition  probability,  we  need  to  find  l/WI2  =  IM0I2  +  +  M0Mf  +  I/W-, I2.  The 

term  IA^,  I2  oc  (eAo/m)2((o/m)2  =  10-13  and  can  be  safely  ignored  in  comparison  with 
IM0I2.  Also,  the  cross-terms  like  M$M:  «  (eAo/m)(a>/m)  =  10-6,  so  that  ignoring  these 
terms  when  calculating  \M\2  will  introduce  an  error  of  one  part  in  106.  Therefore  we 
will  use  I  Ml2  =  IM0I2.  Carrying  out  the  trivial  integration  in  (124)  we  find 

\M,\2  =  l5<4>(/-p+p'-2(Q-P)Ac)!2  (126) 

x  lo^^OycV^)!2  G02  (Q.a,P) 

for  the  specific  harmonic  number  q.  We  note  that  in  order  for  the  energy-momentum 
conservation  relationship  inside  the  6-function  to  be  true,  q  must  be  a  positive  integer. 
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The  total  transition  probability  is  then  'Z\Mq\z.  We  write  the  square  of  the  delta 

<7=1 

function  as 


!8*4)(/-p+p'  -2(qf-|3)/r)  I2  =  -^L^\l-p+p' -(q-2m  (127) 

(2ti)4 


where  V  and  T  are  the  interaction  volume  and  time  (Ref.  15). 

The  transition  rate  per  unit  time  \Mq\2/T  into  a  group  of  final  electron  and 
photon  states  is 


dWq  = 


Vd3\p'\ 

Vd3\Ti ' 

T 

(2k)3 

(203  J 

(128) 


where  the  two  differential  factors  are  the  density  of  final  states  for  the  outgoing 
electron  and  emitted  photon  and  integration  over  p'  and  7  is  implied  (Ref.  19). 
Substituting  (126)  into  (128)  we  find 

dWq  =  7^2  \n^\p')Tt{r)u^(p)\2  (129) 

(2k)  t  t 


X  Gq  (q,a,P)  d3lj^'l  d3l7l 


and  we  see  that  the  dependence  on  V  and  T  disappears,  as  it  should.  Performing  the 
integrations  we  are  left  with 

i-i 


dWn 


e2  ft  m  m 

a 

(2n)2  2  E  £■ 

aft 

\n^\p')rt{r)u{s)(p)\2  (130) 


X  Gq  (Q.a.P) 


which  is  the  differential  transition  rate  for  a  photon  emitted  into  the  solid  angle  element 
dL.  The  derivative  term  in  brackets  is  the  result  of  the  171  =  Q  integration  over  a  delta 
function  whose  argument  is  itself  a  function  of  the  photon  emission  frequency  ft  (Ref. 
19).  We  note  that  (130)  is  subject  to  the  energy-momentum  conservation  relations 
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which  implies 


and  because 


Pu'  =  Pn  +  2(9  -  P)^  -  ^ 

(131) 

P'  =  P  +  2(q-0-T 

(132) 

E  =  E  +  2(9  -  P)co  -  £2 

(133) 

E  =  (\p'\2  +  m2)v*  , 

(134) 

E  itself  depends  on  £2  through  (132).  Equation  (131)  is  Lorentz  invariant;  it  holds  in 
all  inertial  reference  frames. 


The  remainder  of  the  calculation  will  be  easiest  if  we  use  the  "laboratory  frame" 
defined,  for  historical  reasons,  as  that  in  which  the  incoming  electron  is  at  rest.  This 
implies  that  p  =  0  and  p^  =  (0,  im).  Thus  (132)  becomes 

p'  =  2(q-0-T  (135) 


and  the  outgoing  electron  energy  is,  substituting  (135)  into  (134), 


E  = 


4 (q  -  P)2co2  +  £22  -  4(q  -  P)ot)£2cos0  +  m2 


v? 


(136) 


where  0  is  the  angle  between  /r  and  T.  Since  we  know  J<  =  (2 n/X)i,  0  measures  the 
angle  of  photon  emission  from  the  z-axis,  as  shown  in  Figure  30. 

We  can  find  the  frequency  £2  of  the  emitted  photon  by  squaring  both  sides  of 
(132)  to  obtain  the  invariant  equation 


2(9  -  P)p-k  -  I  p  -  2(q  -  P)/  k  =  0  (137) 

where  we  have  used  p2  -  p'2  =  0.  Computing  the  scalar  products  in  (137)  in  the  lab 
frame  and  solving  for  £2  we  find 


_ 2(q  -  P)<» 

1  +  2(9  -  p)— (1  -  cos0) 
m 


(138) 


This  frequency  is  corrected  in  that  it  is  slightly  decreased  from  the  plane  wave 
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approximation,  since  p  «  A§  (1  +  M(kz0)).  The  frequency  is  shifted  by  the  factor 
(1  +  Mkz0),  which  is  the  shift  in  frequency  of  the  incident  laser  light  as  seen  by  the 
incoming  electron,  as  discussed  in  section  A. 


Figure  30.  Coordinate  system  showing  kinematics  of 
Compton  scattering  (after  Ref.  1 5). 


To  compute  a,  (115),  and  p,  (116),  we  compute  the  scalar  products: 

p-k  =  -m o)  ,  (139) 

p'  -k  =  p’-l<  -  E  to  =  co[2(q  -  P)cd  -  Qcos0]  ,  (140) 

ap=  0  ,  (141) 

and 


ap’  =  Xp’  =  -(osin9cos<)>  . 
Substituting  these  into  (115)  and  (116)  we  arrive  at 


(142) 


a  =  eA0 


1  + 


1 


kZr 


sinGcoso 


Q(1  -  cose)  -  m 


(143) 


where  o  is  the  angle  between  the  projection  of  /  onto  the  x-y  plane  and  the  x-axis  (see 
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Figure  30)  and 


e%2 

8com 


1 

0(1  -  cosG) 

kz0 

m  -  0(1  -  cosO) 

(144) 


and  we  see  the  dependence  of  p  on  the  emitted  photon  frequency  O. 

We  are  now  in  position  to  calculate  the  derivative  in  (130).  Using  (136),  (138), 
(144)  and  a  lot  of  algebra  we  find  that 


_a_ 

ao 


(O  +  E'  +  2pco)  = 


2qcom 

EO 


(145) 


We  recognize  this  term,  apart  from  the  harmonic  number  q,  from  ordinary  Compton 
scattering  (see  Ref.  15,  page  228). 

We  can  now  write  (130)  in  the  lab  frame  as 

\a{s,)(p')ye{r)u{s)0)\2  G§(q, a,p)  (146) 

where  a  =  1/137  is  the  fine  structure  constant  (Ref.  15).  The  last  thing  we  need  to  do 
is  calculate  the  term  \ny-eu\z.  Since  the  spin  of  the  initial  and  final  electrons  is 
unobserved,  we  average  over  the  initial  spins  and  sum  over  the  final  spins  with  the 
result 


dW, 


dL 


q  O 
4 n 


O 


q  to 


1  £  lD^><j»')T.£l'>u(s>p)|2  -  E-lJIL  (147) 

2  s.s”  2m  v  ' 


where  we  have  used  "Casimir’s  trick”  in  order  to  compute  the  spin  sum  using  the  trace 
method  (refer  to  pages  191-192  in  Ref.  15). 


Our  final  result,  for  the  emission  rate  of  a  photon  of  frequency  into  the  solid 
angle  element  cfl  around  0,  including  1)  and  c,  is 

dW„ 


dl 


—  =  a— 


4n 


O 

E’  -  me2 

q  to 

2  me2 

^  / 

G5(g,a,p) 


(148) 


Recall  that  for  the  entire  spectrum,  (148)  must  be  summed  over  the  harmonic  number 


96 


q.  We  summarize  our  other  results,  including  D  and  c,  as  well: 


a  = 


2(3  ~  P)<o 


1  +  2(3  -  p)-^(1  -  cose) 
mcd 


E  =  |^4(3  -  +  if Q2  -  4tf(q  -  P)coQcos0  +  m2c4 


Vi 


ti  =  eA0 


1  + 


kZo 


sinQcostj) 


e2An 


P  = 


4mc2 


1  + 


ftQ(1  -  cos0)  -  me2 
1 


kzr 


(1  -  COS0) 


me2  + 


ezA% 


4  me2 


1  + 


1 

kZn 


(1  -  COS0) 


(149) 

(150) 

(151) 


(152) 


and  we  have  used  (138)  to  write  p  in  terms  of  the  scattering  angle  and  laser  frequency 
a>.  Our  Gaussian  mode  correction  (1  +  Vkz0)  is  evident  in  each  of  these  equations 
either  explicitly  or  through  p.  We  see  that  both  the  angular  distribution  (through  G0  in 
(148))  and  frequency  (149)  can  be  altered  by  merely  focussing  the  laser  beam. 

If  we  specialize  now  to  backscattering,  0  =  it,  we  see  that  r|  -»  0  and  the 
Bessel  function  sum  collapses: 


G0(3,T1=0,P)  =  lJZ{n+q)(0)Jnm  =  H)%(P) 


(153) 


since  only  J0(0)  >  0  implies  n  =-q  and  we  have  used  the  identity  (Ref.  20) 

J-nW  =  H  )nJ„(x)  .  (154) 

Therefore  we  find  that  the  transition  rate  becomes 

_  _2 1 

4?(P)  •  (155) 


dWq 

=  a-° 

Q 

E'  -  me2 

dl 

4ji 

q  co 

w  d 

2  me2  , 

The  argument  of  the  Bessel  function  becomes 
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E.  CONCLUSIONS 

In  1935  D.  M.  Volkov  derived  the  Dirac  wavefunctions  for  an  electron  in  an 
electromagnetic  plane  wave.  We  made  a  correction  to  those  wavefunctions  by  solving 
the  Dirac  equation  in  the  external  field  of  a  laser  with  a  realistic,  finite  Rayleigh  length 
z0.  We  then  used  those  wavefunctions  to  find  the  rate  and  angular  distribution  of 
photon  emission  in  this  field.  We  found  that  both  the  angular  distribution  and 
frequency  of  the  emitted  photons  can  be  altered  by  changing  the  Rayleigh  length,  that 
is,  focussing  the  laser. 

Further  research  should  add  the  finite  length  of  the  laser  pulse  to  the 
calculation.  This  will  give  the  emission  lines  a  width  in  frequency  corresponding, 
through  the  Fourier  transform,  to  the  pulse  length. 
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